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•^ Abstract 

■^ We present a bijcctive algorithm with which an arbitrary permutation decomposes 

^_^ canonically into elementary blocks which we call families, which are sets with a specified 

/^ number of ascents and descents. We show that families, arranged in an arbitrary order in 

r \ a sequence, are in bijection with permutations. 

^*^ The permutation decomposes canonically, by inserting parentheses, into a tree having 

i-Q as nodes a class of permutations which we call primitive. Primitive permutations can be 

5^ assembled from very simple data. 

C The data for the trees into which a permutation decomposes can be written in a form 
I— —I similar to the decimal classification of a library. We axiomatize that data. It has a 
structure very different from the permutation which it encodes, with shuffles and pair- 
s' ings instead of reordcrings. These structures arc similar to the fundamental processes in 
j«^ quantum field theory. 

\^ Our main bijcctive structure algorithm gives explicit, additive multinomial formulae 

04 for the number of permutations with given sets of elements under and over the diagonal, 

^~~J or with given ascent and descent values. 

^^ The multinomial expressions obtained this way give a new class of bicolored set statis- 

^^ tics, between set partitions and set compositions, called shifted multinomials. These 

^^ provide for the first time additive multinomial expressions for Eulerian numbers and de- 

. . rangements, as part of a sequence of new combinatorial objects. 

^ These multinomial expressions satisfy inductive relations involving only immediate 
neighbors, similar to the relations satisfied by the Eulerian numbers. 
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1 Preintroduction 

Much of this paper is based on the structure of a family. There are different equivalent ways 
to look at it. First, a family S is a set, together with a specified number k ascent values, or reds, 
and / descent values, or blues, with k,l > 1, and k + 1 = \S\. The information in the bicoloring 
of S can equivalently be given by making 5 into a sequence (ri, r2 . . . , r^, 6;, 6;_i, . . . , 6i), where 
bi < b2 < ■ ■ ■ bi < ri < r2 < ■ . . < Vk, and {ri, . . . , r^} are the reds, the elements colored red 
in the set S, and {bi, . . . ,bi} are the blues, the elements colored blue. We put both kinds of 
information together by viewing 5 as a colored sequence. 



Families 



Permutations are built up by inserting the new (Wghest) element n in a line permutation of n-1, 
in an arbitrary position. Here, 11 is inserted. The element following 11 is called ll's friend. 
The main structure is called a family. Below, 11 tries to fit into the family of its friend. 
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bicolored set. 
At least 
one red and 
one blue; 
All reds 
larger than 
all blues 



A sequence 
family 



The start of a new family 

11 tries, but does not fit in front of its friend 9: 
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(6,8,9,3,2) 

The same family as a sequence. 
The relative position in it 
determines the colors 
reds = ascending values 
(first or larger than previous) 
blues = descending values 
(smaller than previous) 
Reds located first. 
All reds larger than all blues. 




new 
family 



o^ 



oo 



9 would be a 
descent above . 

the ascent 6 3* 




v© 6 



^ remainmg 
fO old family 



(6,8,(11),9,3,2) (6,8,(11,9),3,2) 



Exiting a family 
the highest (new) 
does not fit into the 
family of its friend 
(friend = element 
located to its right.) 



A new family is formed 
by the new element, which 
pulled its friend out. 
A new branch of the 
family tree, parenthesized, 
is started by them. 



The growth of an existing family 
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Entering the family The highest (new) element 
of a blue friend entered the family 

from above. As a red of its blue friend 
between reds and blues (located to its right.) 
(always possible) This creates an extra red. 



of a red friend 
from above 
at the beginning 
(when there is a 
unique red.) 

The only friend positions for which a family can grow are the highest blue and the unique red 



Entering the family The highest (new) 
element enters 
the family of its red friend 
(located to its right.) 
The friend changes color. 
This creates an extra blue. 



The main function of a family S is dynamical. We shall build up permutations of n = 
{1, . . . ,n} in the manner of Euler, by inserting the new element n into the line form of a 
permutation of n — 1, on one of the n positions available. Unless the new n is at the end, 
where it would become a singlet, n is followed in the line permutation by an element m, which 
we call n's friend. The friend m is in a family S, and we try to see whether n fits into that 
family, positioned immediately before m. 

Either S with n inserted in front of m is a family, in which case we say that n fits in S, or 
else n pulls out its friend m, and they start a new family (n, m). If n is inserted at the end 
of S, it forms a new singlet family (n). 

This is the main dynamics of the inner structure of a permutation. The insertion of new 
elements is repeated, as the permutation is built, by inserting new, largest elements, into its 
line description. The tree of possibilities in this process, enlarging families and shrinking them 
as they divide to create new families, is extremely complex, and gives rise to several tree and 
shuffle type structures. Several algorithmic methods, which are the main results of this paper, 
show that the dynamics of family growth and division can be controlled by simple data, or 
reversed to yield bijective structure results on permutations. 

The processes involved are similar to creation and annihilation operators in particle 
physics, or cell growth processes in biology. 

In the paper which follows, we shall study the inner structure of a single permutation. 
Composing two permutations gives very rich structures, not studied here. 

2 Introduction 

Permutations are central to the study of symmetry in mathematics and physics. There are 
many interesting algebraic structures on the linear combinations of permutations. The most 
remarkable of these is the block decomposition of their convolution algebra, namely the study 
of their representations. This has been an active subject of research for over a century, and 
gave birth to modern representation theory. 

In contrast to the above, we shall concentrate in what follows on a new internal structure 
of each individual permutation. Leonhard Euler found the recurrence relations on the number 
of permutations with a given count of ascents and descents, numbers which are now called 
the Eulerian. These have found applications as numbers and operators in many areas of 
mathematics and physics. After three centuries of study of permutations with given numbers 
a of ascents and b of descents, there appear to be no simple additive combinatorial formulae 
in terms of a and b. Such an expression would correspond to a construction of a permutation 
from a list of subsets involving a and b. 

The central result of our study is such a formula. The theory behind it is a way to de- 
compose an arbitrary permutation of {1, . . . , n} into an ordered sequence of mutually disjoint 
subsets Si L) . . . L) Sn in which for every Si, which we call a family, we specify a number 
Oj S {1, . . . , |S'j| — 1} of ascent values, with the remaining number bi = \Si\ — ai counting 
descent values. The algorithm is bijective and provides an explicit way of building the permu- 
tation out of the list (^i, . . . , Sk) which we call a registry of families. The numbers not used 
in the registry become singlet blocks. A transformation, called the cycle transform, sends a 
permutation with given ascent values, descent values and singlets into a permutation with the 
same values which are over, under, and respectively on the diagonal. Our algorithm gives thus 



an explicit construction of the permutations in which such sets, not only their cardinality, are 
specified. 

We show that any permutation has a canonical forest structure, in which nodes of trees 
consist of a special class of permutations, which we call primitive. A primitive permutation 
consists of an active family, with ascent and descent values as described before, together 
with insertion points, called buds, for higher branches. For the tree structure of a permuta- 
tion, primitive permutations play a role analogous to prime numbers with the commutative 
multiplication replaced by a tree structure. 

We give a bijective algorithm for constructing directly the primitive permutations from 
simple combinatorial data such as a sequence of positive and negative numbers - or equiva- 
lently a Dyck path with decorated descents - together with the desired number of ascents of 
the active family. The explicit nature of the data leads to a differential equation with which 
we compute the generating function for the number of primitive permutations with an active 
family of size m and with k buds. 

The data for primitive permutations, which are the nodes of the tree structure of a per- 
mutation, is assembled in a form similar to the decimal classification of a library, which we 
call a decimal code of a permutation. The decimal code, which we axiomatize and show to 
be equivalent to the permutation data, replaces the change of order in the permutation with 
creation - annihilation pairing data and shuffling. 

The description of permutations in terms of sets with a given number of ascents and 
descents gives rise to a study of partition of sets into bicolored subsets. While such partitions 
were first studied half a century ago, we obtain for the first time Eulerian type inductive 
relations for these, involving only neighboring terms, as well as a new type of expressions for 
their count, as sums of shifted multinomials. Such expressions interpolate between bicolored 
partitions and compositions. 

The decomposition of a permutation into a forest is constructed inductively by inserting 
the new highest element in the line description of a permutation, in the manner of Euler. The 
growth of the forest in this process has similarities with biological growth patterns. Small 
nodes, or families, mostly grow while larger nodes tend to shrink in size due to a process 
of reproductive division. This process accounts for the extraordinary complexity of the end 
result of the forest and node structure as well as for the connections to particle creation 
and annihilation of particles in quantum field theory, as well as to possible connections with 
mathematical biology. 

2.1 The library metaphor 

We are in a big library, in which books are labeled by their date of entry. A new book is put 
in a random position every day by a careless librarian assistant. 

The books are arranged on shelves, some, on the left, leaning left, and some, on the right, 
leaning right, at least one on each side. On both sides, numbers should increase toward the 
center of the shelf, and all books on the right must be older than all books on the left. As 
the new, latest book, is put in a random position on a random shelf, there are two cases. The 
new book may fit on the shelf, and this happens in two cases. Either the book is placed in the 
middle, in which case being newest, it will lean left. Or, there is only one book on the left and 
the new book is put against the left wall. Then, the former single left book must lean right. 
In all other cases, the new book does not fit. Then, it picks the book immediately following 



it, and with them a new shelf is started. There are now two ways to proceed, organizing the 
whole library. In the first, the new shelf and the remains of the previous one are arranged in 
the library in a very elaborate way. That allows, with an intricate algorithm, to figure out 
the exact history of the library, just from the positions and contents of the shelves. Moreover, 
every shelf arrangement corresponds to a history. This is our main bijective algorithm. 

Another system is to keep track of the books that were taken away and repositioned, by 
leaving markers for each pair which was taken away, and by assigning these books on their 
new shelf a longer decimal number. The list of these decimal numbers is in mostly ascending 
order, but has books from different shelves shuffled relative to each other, and keeps track 
of the pairings between the books which started a new shelf. This is a decimal code for our 
library, and again it encodes the whole library structure. In this system, which is very far 
from the original permutation, we can think of the pair of books which went missing from a 
shelf, in a grand way, as an annihilation in a universe. As they reappear on their own new 
shelf, they form a new universe. And indeed, the modern type of library, which is the online 
encyclopedia of integer sequences, shows that in the simplest case, when universes disappear 
without creating new ones, on a unique shelf with only two books left on it, the number 
of possible histories matches one to two the number of indecomposable fermionic Feynman 
diagrams with a given number of creation - annihilation pairs. 




Main theorem: 
descent values 



permutations with prescribed ascent and 



3.1 The cycle transform 

Let TT be a permutation of {1, . . . , n}. Extend it with zero at the beginning, i.e., vr(0) = 0. A 
position i is called an ascent (position) if 7r(i — 1) < vr(i), and a descent (position) otherwise. 
If 7r(i — 1) < 7r(i), the value 7r{i) is called an ascent value, else 7r(i) is a descent value. If 
7r{i) > i, then i is called an over diagonal position and 7r{i) an over diagonal value. If 7r(i) < i, 
then i is called an under diagonal position and 7r(i) an under diagonal value. Elements with 
7r(i) = i are called on diagonal. Note that vr maps over/under positions into over/under 
diagonal values. That is not the case with ascent/descent values and positions. Ascents and 
descents are defined between adjacent positions, while the images of ascents and descents are 
not adjacent to each other. 

For the permutation vr, consider the absolute minimum value 1, with 7r(ii) = 1. Let 
mi = 1, and let m2 = minj>jj7r(i), with 7r(i2) = rn2- Continue similarly the sequence of 
successive minima nii < m2 < in^ < . . ., and the corresponding successive minima positions 



ii < 12 < is < ■ ■ ■, with TT{ik) = ruk- We call the intervals {1, . . . , ii}, {ii + 1, . . . ,12}, ■ ■ ■ the 
blocks of the permutation, with ii,i2,--. the block anchors, and mi,m2, ■ ■ ■ the block 
anchor values. A block consisting of a single element will be called a singlet. The similarly 
defined successive maxima, not used here, are called salients in the literature. 

A remarkable transformation was discovered half a century ago by Alfred Renyi, and called 
an "unusual transformation" by Donald Knuth. It is the cycle transform, and is defined as 
follows. 

Decompose a permutation a into cycles, arranged so that the minimum of each cycle is 
on the last position, and so that these minima are in ascending order. Write now the list of 
elements of these cycles, and read it as a line permutation n. For instance, if a, written as 
a line permutation, is 524361, then its cycles ordered as above are (561) (2) (43). The cycle 
transform vr is 561243 in line notation, with successive minima 1,2,3, and blocks 561,2,43, 
which are the cycles of a. a has 4,5,6, respectively 2, respectively 1,3, as over/on/under 
diagonal values, which are also the ascent values, singlet blocks, and, respectively, descent 
values of vr. 

We call the permutation vr the cycle transform of a. Remark that the minima of cycles 
of a become the successive minima of tt, and the cycles of a are the blocks of vr, as defined 
above, a is the inverse cycle transform of vr, and is obtained by interpreting each block 
between successive minima positions of vr as a cycle of a. 

The on diagonal elements of a are the singlet blocks of vr. The over /under diagonal values 
in o are the ascent/descent values of vr, which are not singlets. Depending on the problem 
studied, it is sometimes convenient to view singlets as ascents, while in other applications 
they are set apart in statistics. 

The structure studied in what follows will concern the ascent/descent values and the blocks 
between successive minima of a permutation vr. We shall use the inverse cycle transforms of 
these to over/under diagonal values and cycles only for the interpretation of the results. The 
cycle transform as a three dimensional picture is shown below. 
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3.2 Derangements, Eulerians, and Eulerian lift 

We call a permutation a a derangement if it has no diagonal elements, or cycles of length 
one. We denote by Da,h the number of derangements with a values over diagonal and h values 
under diagonal. The cycle transform sends these derangements into permutations without 
singlet blocks, with a ascent values and h descent values. 

Extend a permutation vr with zeroes to (0, 7r(l), 7r(2), . . . , 7r(n), 0). Euler studied the num- 
ber of permutations with a given number a of ascent adjacent pairs and a number h of descent 
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adjacent pairs in that sequence. The number of such permutations are the Eulerian numbers, 
denoted in the hterature by A{a + b,a), and here by Ea^t, with '^a+b=n+i ^a.,b = nl. 

Given a permutation vr of {1, . . . ,n} with a ascent values, including the singlet blocks, 
and b descent values, as counted by Euler, we define a new permutation vr of {1, . . . , n + 1}, 
which we call the Eulerian lift of vr, by letting Tr{i) = 7r(i) + 1 for i < n and 7f(n + 1) = 1. 
Note that vr has a single block, and the number of its ascent and descent values is the same 
as the number of ascent and descent pairs counted by Euler. For instance, vr = 312 has two 
ascent pairs 0312 and two descent pairs 31, 20 in 03120. Its Eulerian lift is vr = 4231, and has 
in 04231 ascent values 4,3 and descent values 2,1. The inverse cycle transform of tt is the 
cycle (4,2,3,1), written 4312 as a line permutation. It has 4,3 as over diagonal values, and 
1, 2 as under diagonal values. 

The transformation which corresponds to the Eulerian lift by the inverse cycle transform 
is daisy chaining. Let o" be a permutation of {1, . . . , n} and let ci, . . . , Cfc be the list of its 
cycle heads, the smallest elements in each cycle, arranged increasingly. Let a' denote the 
permutation of {0,1,..., n} defined by a'ci = itcj+i for i < k, a'ck = 0, a'O = ac\ and 
^' 2 = <^i otherwise. Now lift {0, . . . , n} to {1, . . . , n + 1} to obtain a. Daisy chaining puts in 
bijection permutations of n with single cycle permutations of n + 1. 

The Eulerian numbers E]^^i count thus the number of permutations of single cycle permu- 
tations with k elements over and I elements under the diagonal. 

3.3 Patterns and shuffles 

In the sequel, we shall make use of patterns and shuffles. For a finite F C M, an injective 
function / : F — )• M has a pattern which is a permutation of (1, . . . , |-F|), obtained by 
composing / on both sides with order preserving maps. 

Given an injective function / : F — )• M on a finite F C M we call / a shuffle of /i, . . . , /„, if 
there is a partition F = \}^Fi so that the restriction of / to any Fi has the same pattern as 

/.. 

More generally, for an injective function / : F = IJ^ F^ — t- M with Fi not necessarily 
mutually disjoint, we call / a shuffle of /i, . . . , /„ along (Fi, . . . , F^ if the restriction of / 
to any Fi has the same pattern as /j. 

3.4 Families and registries 

In what follows, we shall develop structural results which will put a permutation in bijection 
with a composition, an ordered partition, of part of its positions into bicolored sets, called 
families, with the remaining elements singlet blocks. The algorithm will construct explicitly, 
and thus count, permutations with given sets of ascent/descent values, or equivalently, via the 
inverse cycle transform, of given sets of over/under diagonal values or over/under diagonal 
positions. 

Recall that a set family 5 is a subset of {1, . . . , n} with a color, red or blue, assigned to 
each element, with r G {1, . . . , IS"! — 1} reds, and blues, such that any red element is larger 
than any blue element. Thus, a set family is completely determined by its underlying set S 
with > 2 elements, together with the number r € {1, . . . , l^l — 1} of reds. 

A sequence family is a sequence of integers (oi, . . . ,a/c, 6;, . . . , 6i) with A;, / > 1 and 
6i < 62 <•••<&«< tti < ^2 <•••< Ofc- The elements ri, r2, . . . , r^ are called ascent values. 



or reds, and the elements bi,b2, ■ ■ ■ ,bi are called descent values, or blues, and we color the 
sequence as (ai, . . . , a^, 6;, . . . , 61). Note that the reds and the blues in a nonsinglet sequence 
family are completely determined by their relative position in the sequence, and all reds have 
values larger than all blues. Thus the subjacent bicolored set of a sequence family is a set 
family. 

Conversely the set family {61, . . . ,bi,ai, . . . , a^} with 61 < • • • < 6/ < oi • • • < a^ deter- 
mines the sequence family (ai, . . . , a^, 6/, . . . , 61), or (ai, . . . , a^, 6;, . . . , 61). We shall consider 
all these as set and sequence forms of a same family. 

Given a subjacent set S with > 2 elements, there is precisely one family with a given 
number r G {1, . . . ,\S\ — 1} reds. 

We shall also sometimes consider singlets as families (colored black in the pictures). 

A family registry is an ordered sequence of mutually disjoint set families. For various 
purposes, we may skip the singlet families from registries. 

Examples 

S = {2,3,5,6,7} colored with 2 reds is the set family {2,3,5,6,7}, and the sequence 
family (6, 7, 5, 3, 2) colored as (6, 7, 5, 3, 2), or 67532. 

A set with four elements {1234} can be organized as a family in 3 ways: 4321, 3421, 2341. 
These are all the patterns of families with four elements. 

The colored sequences 58431 and 752 and 9 are families while 58341, 58631, 752 and 58 
and 531 are not. 

4 The main counting theorem 

Our main permutation construction result is the following. 

Theorem 4.1. The permutations with a given set A of elements over the diagonal, a given 
set of elements B under the diagonal and the rest on the diagonal are in bijection with the 
number of compositions A = AiU ■ ■ -L) A^, B = BiU ■ ■ -U Bi^ with min Ai > max Bi for all i. 

An example 

Consider the subsets A = {2,4,5} and B = {1,3} of 5 = {1,2,3,4,5}. We want to 
construct all the permutations with A as ascent values and B as descent values, or equivalently 
all the permutations with A as over diagonal values and B as under diagonal values, or 
equivalently all the permutations with B as under diagonal positions and A as under diagonal 
positions. 

Consider the compositions A = Uj^jAj, B = Uj^jBj with min^j > maxBj for all j. 

There are 6 such compositions, written as vectors 

'n ^ o)^!?) ^ {im ^ 0(3) ^ ©^(T) ^ a>© - ©K? 

A pair Aj > Bj forms a set family Aj U Bj with Aj as reds and Bj as blues. The families 
given by the compositions are arranged as registries, i.e. ordered lists of families, as follows. 

((21), (453)), ((453), (21)), ((241), (53)), ((53), (241)), ((251), (43)), ((43), (251)). 

By our reverse algorithm these registries are mapped into 6 permutations vr of 5 with A = 
{2,4,5} as ascent values and B = {1,3} as descent values, respectively (using the block 



algoritni) . 

21453, 24531, 24153, 25341, 24351, 25143, 

The blocks between successive minima of these are respectively 

(21)(453), (24531), (25341), (241)(53), (251)(43), (24351). 

These blocks give the cycles of the inverse cycle transform. Written in line notation, the 
inverse cycle transforms are 

21453, 24153, 25413, 24513, 25431, 24531. 

In these, the elements {2, 4, 5} of A and respectively {1, 3} of B, are the values over, respec- 
tively under, the diagonal, as desired. 

We now replace each permutation a above by a~^. This will make A = {2,4,5} and 
B = {1,3} the sets of positions mapped under (u) respectively over (o) the diagonal. That 
is, we obtain the 6 permutations with the pattern (ououu). 

21534, 31524, 41532, 41523, 51432, 51423 

In our study of the structure of a permutation, we shall decompose it with nested parenthe- 
ses, such that the outermost parentheses are its blocks (between successive minima positions). 
The contents of each parenthesis, in which the included parentheses and their contents are 
removed, form a sequence family. This forest structure, in which each block is a tree, is 
canonical, and is called the parenthesized form of the permutation. 

The parenthesized forms of the permutations vr, which refine their block decompositions 
as written before, are respectively 

(21)(453), (2(453)1), (2(53)41), (241)(53), (251)(43), (2(43)51). 

The above parenthesized forms show that the compositions of the ascent/descent sets A, B 
from which we started, written as sequence families 

{2,4,5,1,3} = (21)U(453) = (453)U(21) = (241)U(53) = (53)U(241) = (251)U(43) = (43)U(251) 

appear, very transparently, in the very body of the corresponding permutations, inside the 
parentheses. 

Our direct algorithm gives, starting from these permutations, an ordering of the families 
which appear in parentheses as a registry of families. This way we recover the compositions 
of A, B from which we started. 

5 Bicolored partitions and their Eulerian type recurrence re- 
lations 

We shall count Eulerian numbers Ea^b^ derangements Da^bj ai^d similar classes of bicolored 
permutations, by a new type of statistics of bicolored sets, which we now describe. 

Let k,l G {1,2,...}. Decompose as vectors {k,l) = (/ci,/i) + ••• + {kr,lr) into an (un- 
ordered) partition, with fcj, li S {1,2,... }. Let ri, ..., r^ be the number of identical parts of 
each type, with r = ri + ■ ■ ■ + rp. 
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The multinomial sum 

>P {k + iy. 






(A;i + /i)!...(A;, + /,)! 



counts the decompositions of the set {1, . . . , A: + /} into r subsets of sizes ki + li, . . . ,kr + Ir- 
The sum 

^H ^ V (^ + 0' ^__ 

^^^ ^{ki + h)\...{kr + lry.n\...rp\ 

counts the corresponding set partitions, i.e. the decompositions in which the order of the 
subsets is neglected. Note that subsets with the same cardinality ki + li = kj + Ij, in which 
the set labels [ki, li), ikj,lj) are different are counted as different in the partition. We call the 
numbers Nj^ j the fundamental r— part multinomials. 
If we now take 

HivM = V ^^^^^^ '-^ 

■ '='' ^(A;i + /i)!...(A:, + /,)!ri!...rp!' 

we count the above as set compositions, or ordered partitions. More generally, we shall count 
shifted partitions, defined by 



ns= e ir+s)iNt}= j: e 



(A:i + /i)!...(/c, + /,)!ri!...rp!- 

r>inax(— s,l) r>max(— s,l) 

Here s S {. . . , — 1, 0, 1, . . .} is a shift, and we take negative factorals to be zero. We call the 
numbers Njfl the shifted multinomials of shift s. 

Our main enumeration result is the following. 

Theorem 5.1. The Eulerian numbers have the multinomial expression 
E..^Nt,^'-^ ^±i^ i^^ 



' ^^ {ki + h)\...{kr + lry.ri\...rpl' 
The derangements have the expression 

n =,v(o) = V (^ + 0' ':L_ 

^'' '^'^ 2^{ki + h)\...{kr. + lr)\ri\...rpV 

summing over the unordered partitions {k, I) = {ki, h) + ■ ■ ■ + {kr, Ir), with ki, U G {1, 2, . . . , }, 
ri, . . . , Tp the number of identical parts of each type, and r = ri + • • • + r^ 

We shall show, with a method which relies on the decomposition of partitions into set 
families, the following Eulerian type recurrence relations, between the fundamental and shifted 
multinomials. 

Theorem 5.2. 

<! = kN^i^^ + lNtl^i + {k + l-l){Nt:^,_,-rNtl,_,) 

nS = kN^:l^ + lNi%, + {k + l-l){NizZ'\ + {^ + mtX^^^^^ 
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Theorem 5.3. For negative shifts s < 0, using compositions (ordered partitions) with r > —s 
parts 

{k,l) = {ki,li) -\ h {kr,lr), 

where fcj, /j G {1, 2, . . .}, letting Ui = ki + k, the shifted multinomials N^^l decompose as 

k,l - Z^y m - 1 )\ 712-1 

n^s-1 H h n^ - 1\ / n_s -\ \- Ur - 1 

This generahzes the expression of multinomials in terms of binomials. The identity de- 
composes the shift s < into |,s| binomial shifts of —1 (and does not work for other such 
decompositions of the shift) . 

5.1 Combinatorial interpretation 

Shifted multinomials with shift count compositions of a bicolored set with a reds and b blues 
into parts Si which have each a number of reds Oj > 1 and blues bi > 1 , \Si\ = Oi + bi. 

For a negative shift — s < 0, s of the parts are on fixed positions. For instance let 5i be 
the part containing the smallest element 1, 5*2 the part containing the smallest element not in 
Si, and so on. Then N^^l counts compositions which keep Si, . . . , Ss, in order, on the first s 
positions. In the case of our algorithm which maps permutations into ordered lists of bicolored 
sets, Si, . . . ,Ss are in order on the last s positions. In particular single block permutations 
are mapped onto compositions with Si last. As we map permutations of n, with the Eulerian 
lift, onto single block permutations of n + 1 and count Eulerian numbers on the image, we 
obtain expressions of Eulerian numbers in terms of shifted multinomials with shift — 1 . 

For positive shifts s > 0, shifted multinomials count compositions in which s dummy parts 
(markers or jokers) are inserted as parts. Thus positively shifted multinomials are divisible 
by s!, accounting for the relative order of the extra parts. 

5.2 Examples 

The above formulae show that the Eulerian number E^^-^ is given by the sum over the partitions 

(3, 3) = (2, 2) + (1, 1) = (2, 1) + (1, 2) = (1, 1) + (1, 1) + (1, 1) 

6!0! 6! 1! 6! 1! 6! 2! 1 

^^'^ = 6! l! + 4!2! Im + M im + 2!2!2! 3! = ^ + ^' + '° + ^° • 3 = ^^- 

The fundamental multinomials are A'^ 3 = 1; -^3 3 ~ ^^' -^3 3 ~ -'-^• 
From these, we compute the Eulerian 

^3 3 = oiN^^j^ + V.N^^l + 2!ivg = 1 + 35 + 2 X 15 = 66. 

From the same fundamental multinomials, with different coefficients, we compute the de- 
rangement number 

1)3 3 = l!iv]^l + 2!^^^ + 3!ivi^J = 1 + 2x35 + 6x15 = 161. 
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The Eulerian number E^^^, also has the following expression, defined in terms of composi- 
tions of (3,3), 

(3, 3) = (2, 2) + (1, 1) = (1, 1) + (2, 2) = (2, 1) + (1, 2) = (1, 2) + (2, 1) = (1, 1) + (1, 1) + (1, 1) 

with sums 

6=2+4=4+2=3+3=3+3=2+2+2 

is, after shifting the size of the first parts by —1, as sums of multinomials 

Eulerians, derangements, multinomials, and their relations are illustrated in the following 
picture. 

SfflFTED MULTINOMIALS, [BICOLORED PARTITIONS], AND THEIR RECURRENCE RELATIONS 

- derangements ^ Eulerian - 

shift +1 shift shift -1 [}] shift -2 shift -3 

2 2 II [partitions] I I 

n part] [1] [1] 

' '" ' ' ' ' [}][!] [}] ^s^A. 

2 62 62 2 I 21 21 I I II II I l^ partsj ^q 10 

[I] [I] [I] [I] [10] [10] 

2 152 572 152 2 I 51 161 51 I I 26 66 26 I 25 50 25 [partitions] 15 

[1] [1] 01 [1] [1] [25] £35] [25] 3 parts] USl 

2 338 3068 3068 338 2 I 113 813 813 113 I 1 57 302 302 57 I 56 196 196 56 ^ i* ^ 105 105 

[1] [1] [1] [1] [1] [1] [56] [91] [91] [56] [I05][I05] 

572 = 3-62+3-62+2-5-20 161 = 3-2I+3-2I+I-5-7 66 = 3II+3II+0-5-4 SO = 3I0+3I0-I-5-3+5I 

1 = 3- 1+3- I-I-5I 35 = 3I0+3I0-2-5-3+5I 15 = +5-3 

572 = 2!I+3!-35+4!I5 161 = I!I+2!-35+3!I5 66 = 0!I+I!-35+2!I5 50 = 0!-35+I!I5 I5 = 0!I5 

5.3 Comments 

In the middle of the last century, colored partitions were studied by analogy with the rich 
structure of number partitions, but the interest waned without the realization that bicolored 
partitions were fundamentally related to permutations. In particular our shifted multinomials 
appear to be a new concept, missing at that time. In our case, the shifts in factorials are due 
to linear relations in the Quantum Field Theory model which motivated the combinatorics. 
The shifted multinomials Nj^l satisfy Nj^^ = Ej^^i, the Eulerian numbers, and Nj^^ = 

D^i, the derangement numbers. A'^^ ^ appeared in the count of certain hyperdeterminants. 

(— fc— 1) 
The boundary cases Nj^ ^ are associated Stirling numbers of the second kind. 

(s) \r] 

Except for these few, according to the oeis database, no other numbers Nj^ ( , Nj^ j which 
we described above appear to have been studied before. 

The connection between the Eulerian Nj^ ^ and derangement N^, I numbers, both obtained 

from the same fundamental numbers N^ j with different factorial multiplicities appears to be 
new as well. 

In the literature, there are relations for bicolored partitions involving an arbitrarily large 
number of terms. Our Eulerian type relations, which involve only immediate neighbors, 
appear to be new. 
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6 Recurrence relations for multinomials. Algorithms and proofs. 

We now prove the Eulerian type relations 

<1 = kNtl, + lNtl,, + {k + l-l){Nt,%-rNtl,,^,) 
combinatorially, by providing a bijection between the two members. The equality involves the 

\r] 

unordered partitions number N^, j, which counts bicolored subsets which are mutually disjoint, 
each with at least one element of each color, with a total of k respectively / elements of each 
color, and having the elements labeled by the numbers 1,2, . . . ,k + 1. The sequence of sets is 
considered unordered, and inside each set the assignment of colors to elements is ignored. In 
other words, we look at an unordered partition of 

{1,2,...,A; + /} = S'iU---U5',., 

and numbers ki,li, . . . ,kr,lr S {1)2,...,}, with ki + k = \Si\ and ^^ ki = k, ^^ k = I. We 
color each set Si as a set family, the highest ki numbers in each Si red and the other li blue. 

r\r] \r] 

In what follows we shall denote by Mj^l the set, having cardinality Nf^i, of all partitions of 

[21 

{1,2,...,} into nonsinglet set families with a total of k reds and / blues. For instance, Afj^ 2 
consists of the partitions 

{1,2}U{3,4},{1,3}U{2,4},{1,4}U{2,3}, 

or more simply 

12U34, 13U24,14U23. 

We write the induction relation in a positive form, as 

We shall describe a bijection between sets counted by the two members above. A'^^ ^ counts 
partitions of {1, . . . ,k + 1} into r families, with a total of k reds and I blues. For such a 
partition V we consider 3 cases, depending on the family S of its largest element k + l, which 
being largest is colored red, 

Case (A) The family S is small, i.e. it has two elements. 

Case (B) The family 5 has > 2 red elements (including k + V). 

Case (C) The family S has A; + / as the unique red element, and > 2 blue elements. 

The second term r{k + / — l)Nl_^f^_-^ counts triples {V,m,S) where 7^ is a partition of 
{1, . . . ,k + l — 2} into r families, with a total of A; — 1 reds and l — l blues, m & {1, . . . ,k + l — l}, 
and 5" is a choice of one of the r families of V. 

We consider two cases. Decompose S = RU B into its red and blue elements, where r > b 
for any r & R and b £ B. Now either 

Case (D) m < raaxB or 

Case (E) m > maxS. 

The corresponding cases counted by the right side member are the following. 

Case (A') The first term [k+l — l)N^_-^ ^_^ counts pairs {V, m) between a partition V into 
r — 1 families having a total of A; — 1 reds and I — 1 blues and a number m G {1, . . . , A; + / — 1}. 
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The second term kN^\_^ counts pairs ("P, m) between a partition V of {1, . . . , A; + / — 1} 
into r families having a total of k reds and I — 1 blues and a number m which is among the k 
reds in V. Let S be the family containing m. We distinguish two cases. 

Case (C) The family S has m as its unique red element, or 

Case (E') The family S has > 2 red elements, including m. 

The third term lNj^_^ i counts pairs {V, m) between a partition V oi {1, . . . , k + 1 — 1} into 
r families having a total of k — 1 reds and I blues and a number ?ti which is among the I blues 
in V. Let S" be the family containing m. We distinguish two cases. 

Case (B') m is the largest blue in S, or 

Case (D') m is not the largest blue in S. 

We now describe a bijection between each case A, . . . ,E and the corresponding case 
primed. 

Case (A) 7^ is a partition of {1, . . . , A; + /} into r families, with a total of k reds and I 
blues. Its largest element A; + Hs in a small family S = {m, k + I}. Remove the family S from 
the partition V. Lower all labels > tti by 1, to obtain a partition V' of {1, . . . ,k + I — 2} into 
r — 1 families having a total of A; — 1 reds and / — 1 blues. The pair (V', m) is in the Case (A'). 

Case (A') Given a pair (P, m) between a number m € {1, . . . , A; + / — 1} and a partition V 
of {1, . . . , A; + Z — 2} into r — 1 families having a total of A; — 1 reds and I — 1 blues, increase all 
labels > ?n in 7^ by 1 and append to V the family {ttt.. A; + /}. We obtain this way a partition 
V' of {1, . . . , A; + /} into r families having a total of k reds and / blues, with its largest element 
A: + / in a small family, which is the Case (A). 

Case (B) We have a partition V of {1, . . . , A; + ?} into r families having a total of k reds 
and I blues, with the family S of the largest element k + I having > 2 red elements. Let 
m denote the largest blue in S. Remove the element k + I. The remaining partition V' of 
{1, . . . , A; + / — 1} into r families has a total of A; — 1 reds and I blues. The pair (7^', m) is in 
the Case (B'). 

Case (B') We have a pair {V\ m) between a partition V of {1, . . . , A; + / — 1} into r families 
having a total of A; — 1 reds and I blues, and a number m marked blue, which is the largest 
blue in its family S. Insert the element k + I colored red into S to obtain a partition V of 
{1, . . . , A; + /} into r families with a total of k reds and / blues, with the family S of its largest 
element k + I having > 2 red elements. Then V is in the Case (B). 

Case (C) We have a partition V of {1, . . . ,k + 1} into r families having a total of k 
reds and / blues, with the family S of the largest element k + I having no other red element. 
Remove the element k + I. Let m be the remaining largest element in S. Change the color 
of m from blue to red. The remaining partition V' of {1, . . . , A; + Z — 1} into r families has a 
total of k reds and / — 1 blues. The pair ("P', m) has m as the unique red in its family, so is 
in the case Case (C). 

Case (C) We have a pair {V\m) between a partition "P' of {1, . . . , A: + / — 1} into r 
families having a total of k reds and I — 1 blues, and m which is a number in {1, . . . , A; + / — 1} 
marked red in V' , with m being the unique red in its family S. Insert k + 1 into S and change 
the color of m from red to blue. The corresponding partition P of {1, . . . , A; + /} into r families 
has a total of k reds and I blues. The largest element A: + / is the unique red in its family. 
The partition V is in the Case (C). 

Case (D) We have a triple {V,m,S) where l<?n,<A; + ^ — 1, Pisa partition of 
{1, . . . , A; + / — 2} into r families, with a total of A; — 1 reds and / — 1 blues and S" is a family 
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of V. The family S = R^J B is decomposed into its red and blue elements, and m < maxi?. 
Raise the labels of elements > ?7i in "P by 1 and insert m colored blue into S. We obtain a 
partition V' of {1, . . . ,k + 1 — 1} into r families, with a total of A; — 1 reds and I blues. In it 
m is not the largest blue in its family. The pair (V',m) is in the Case (D'). 

Case (D') We have a pair {V' , m) between a partition V' of {1, . . . ,k + I — 1} into r 
families having a total of k — 1 reds and I blues, and m which is a number in {1, ... ,k-\-l — 1} 
marked blue in V' which is not the largest blue in its family S' . Remove m and lower the 
elements > mhy 1 to get a set 5 and a partition V of {1, . . . ,k + l — 2} into r families having 
a total of A; — 1 reds and / — 1 blues. Then the triple (V, m, S) is in the Case (D) 

Case (E) We have a triple {V,m,S) where l<m<k + l — 1, Visa partition of 
{1, . . . ,k + I — 2} into r families, with a total of A; — 1 reds and / — 1 blues and S is a family 
of V. The family S = RU B is decomposed into its red and blue elements, and m > maxB. 
Raise the labels of elements > m in P by 1 and insert m colored red into S. We obtain a 
partition V' of {1, ... ,k + l — 1} into r families, with a total of k reds and I — 1 blues, so that 
m is not the unique red in its family. The pair {V' , m) is in the Case (E'). 

Case (E') We have a pair {V\ m) between a partition V' of {1, . . . , A; + / — 1} into r 
families having a total of k reds and I — 1 blues, and m which is a number in {1, . . . , A: + / — 1} 
marked red in V\ so that m is not the unique red in its family S' . Remove ttt,, lower the 
elements > ?n, by 1 to get a family S and a partition V of {1, . . . , A; + / — 2} into r families 
having a total of A; — 1 reds and / — 1 blues. Then the triple {V, m, S) is in the Case (E) 

This ends the proof. 

We choose an arbitrary integer s and prove the recurrence relation for s-shifted multino- 
mials. We multiply the relation obtained 

^tl = kNtl, + lNtl,, + {k + l-l){Nt,%,-rNtl,,,). 



with (r + s)!, where negative factorials are taken to be 0. We obtain 

(r + s)!ivg„^ + /(r + s)\NJ[l^^i + 
+{k + l- l)((r + s)!ivtili - <"■ + ^y-^t-u-i) 



(r + s)!ivM = k{r + sy.NPi^^+l{r + s)\NJ[l^i + 



which we write as 

+(k + l- l)(((r -l) + is + l))!ivt"ili - (^ + (^ + mr + s)lN^^l,^,_^) + 
+(A: + / -!)((. + l)(r + .)!iv]rii,_i). 



We now sum over all r > max(— s, 1). The sum on the second row is telescopic with respect 
to r, and we are left with 

(A; + / - l)(max(-., 1) + s)!ivS;^/_7''^''l = {k + I - l)(max(. + 1, 0))!iv|;!^;/„7-''°^' 



If s > -1 thenivi;^^^/_7-''°^^ = ^k-u-i = 0- Else if s < -2 then (max(s+l, 0))!ivS;^/_7-''°^^ 
ob 

is) 



^k-i i-i- ^^ obtain in all cases 



KJ = kKJ_, + lNi:!,^, + {k + l-l){Njr-\ + {s + l)Nl\,_,) 
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Note that when s > — 1 the last relation simplifies to the positive 

Njfl = kNlfl, + lNi%^, + (A: + / - 1)(. + l)iVa,,_i, 
In particular this gives for s = —1 the relations satisfied by Eulerian permutations 

and for s = the relations satisfied by derangements. 

7 The structure of a permutation: families as building blocks 

7.1 The entry /exit lemma 

Recall that a nonsinglet family S written as a bicolored set is S = {bi,b2, ■ ■ ■ ,bi,ai,a2, ■ ■ ■ ,ak}, 
with > 1 reds and > 1 blues, in which any red is larger than any blue. The same family written 
as a sequence is S = {ai, . . . ,ak,bi, . . . ,bi) with k,l > 1 and 6i < 62 < • • • < ^/ < ^i < ^2 < 
• • • < flfe, or as a colored sequence S = (ai, . . . , a^., bi, . . . , 61). 

Definition. Let S be a nonsinglet family, written as a sequence. Let ?7i G M be such that 
m > max 5. Denote by 5" the sequence obtained by inserting m in the sequence S on any 
position other than the last at the end of S. Let x £ S he the element located immediately 
after m in the sequence S' , which we call the friend of m in S. We call x an exit element 
of 5 (or say that x is in an exit position in S) if 5' is not a family. 

Remarlcs. A family with 2 elements, which we call a small family, has no exit elements, 
while a family with > 3 elements, which we call a large family, remains a family after the 
removal of any of its elements. 

In the algorithms which follow, the friend 2; of m as above will be pulled out of its family 
S and will start a new sequence family, {m,x). From the perspective of the old family S, 
[m, x) will be called a bud, as it starts a new branch of the family tree. From the point of 
view of the newly started family, [m,x) will be called the family root, with x called the 
family founder and m the family cofounder. 

Definition Let 5 be a nonsinglet family, written as a bicolored set. Let 2; G N be such 
that X ^ S. Let S' = SU {x}, bicolored so that it is a family and the restriction of its coloring 
to S is the coloring of S. We say that x had a regular entry in S\ or was regularly 
inserted into 5', if the color of x in S' is determined as follows. In all cases except for 
one, the color of is determined by the requirement that after x enters, 5 is a family. The 
only exception is when the value x is between the highest blue and the lowest red. In that 
case, placing x at the beginning, as a red, or between the reds and blues, as a blue, will 
both give families. The algorithms in which this process is used require that after x enters 
5', then x, as a friend should be possible to extract by a higher element placed before it in 
S. This is possible only in the case in which between these two choices x enters as a red in 
front, and we shall call this case a regular entry. Thus a regular entry is characterized by the 
fact that it can be reversed by an exit, which is the content of the entry /exit lemma. Let 
S = {61, 621 • • • 5 h, ai, 02, . . . jflfe}, with 61 < 62 < • • • < ^i < fli < 02 < • • • < fflfc as above. If 
X <bi then x is colored blue in S' . If x > 6/ then x is colored red in S' . 
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Remark If 2; < bi, respectively x > ai then x must be colored blue, respectively red in 
order for S' to be a family. The only multiple choice of color is when 6; < x < ai and in 
that case a regular entry will color x red, in view of the result which follows, which connects 
the horizontal (sequence order) of a family in sequence form with the vertical structure of the 
same family as a set family. 

Theorem 7.1. (the entry /exit lemma) 

If S is a nonsinglet family, x £ S is in an exit position, and S' denotes the family S with 
X removed, then the regular insertion of x into S' gives back S. 

If S is a nonsinglet family, and x ^ S is inserted regularly into S to obtain a family S' , 
then X is in an exit position in S' . 

Proof. The only nonexit positions in a family S are the highest blue and the unique red. 

If 5 = (ai, . . . , Ofc, bi, . . . , bi) with bi < b2 < ■ • • < bi < ai < a2 < ■ • • < a^ and m > max 5, 
then <S" = (ai , . . . , Ofc , m, 6; , . . . , 5i ) with bi < b2 <•••< bi < ai < 02 <•••< Ok < m is a 
family, in which the friend of ?7i is 61, the highest blue. 

In case k = 1, i.e. if there is a unique red, then S' = (m,ai,bi, ... ,61) with 61 < &2 < 
• • • < 6/ < ai < ?n is a family in which the unique red ai in S has changed color to blue in 
S", and m is now the unique red in S' . It is easy to check, and left to the reader (see also the 
illustration) that if the friend of m is any element other than the highest blue or the unique 
red, then S' is not a family, so that a position in S is an exit position iff it is not the highest 
blue or the unique red. Note that in a small family l^l =2, the 2 elements are the unique red 
and the highest blue, so there are no exit positions. 

We now check that if a; € (S* is in an exit position and is taken out of S to obtain S' , then if 
X is reintroduced regularly in S' then x recovers the color it had in S, and thus x recovers its 
position in S as a sequence family. If x was blue, but not the highest blue, it will enter below 
the highest blue and be colored blue. If x is red, but not the unique red, then the highest 
blue in S will remain the highest blue in S", so x will enter above the highest blue in S' and 
with a regular entry will be colored red as it was before. 

For the second part of the statement, remark that in a regular entry of x in S, either x is 
below the highest blue of S and will be colored blue, but will not become the highest blue, or 
X is above the highest blue of S and will be colored red. In the latter case x will not become 
the unique red, since S had at least one red already. Thus in both cases x will be in an exit 
position in S' = SU {x}, since the only nonexit positions are the highest blue and the unique 
red. D 

Recall that a nonsinglet family F written as a bicolored set is {bi,b2, ■ ■ ■ ,bi,ai,a2, ■ ■ ■ ,ak}, 
written as a sequence is (ai, . . . , Ofc, 6;, • • • , &i) or as a colored sequence S = {ai, . . . , Ok, bi, . . . , 61), 
with k,l > 1 and bi < b2 < ■ ■ • < bi < ai < a2 < ■ • • < a^; its subjacent set is 
{61, . . . ,bi,ai, . . . , Ok}- A singlet family is as a set {s} and as a sequence (s). 

Definition A (family) registry 7^ is an ordered list, or sequence, of mutually disjoint 
families TZ = {Fi, . . . , Fk). The set S = [j^Fi is called the range of 7^. Let (p be the order 
preserving bijection 5— >•{!,... jlS'j}. The pattern TZ' = {F[, . . . , F^) of 7?. is the registry in 
which F- = (j){Fi) and (j) preserves colors. 

A stripped registry is a registry with no singlets. From a full registry one obtains a 
stripped registry by removing the singlets. 
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A registry TZ is decomposed into blocks by successive minima in the same way as a per- 
mutation. Namely, let mi be the minimum of all members of all families in TZ, found in a 
family Fj^ Let m2 be the minimum of all members of all families located after Fi-^ , found in 
a family Fi^, etc. With zq = we have now registries 7^^'^ = (Fj;_^+i, . . . ,FjJ called block 
subregistries which concatenate to TZ. For a block of a permutation, its last element, which is 
also its minimum, is called the block anchor. Similarly, for a block subregistry, the minimum 
of all the families in it is called the anchor. 

Suppose that 7^ is a registry with support n in which the only singlets are singlet blocks. 
Let TZ' be the registry from which all the singlet blocks have been stripped, i.e. removed. Let 
S denote the complement of the range of TZ' in n. Insert the elements of S as singlets into 
TZ' successively, starting with the smallest one, in increasing order. If the smallest singlet is 
(1), insert it on the first place. Else compute the blocks of TZ' , which end in the sequence 
of successive minima. Insert the new singlet (s) on the position immediately to the right of 
the last block which has minimum > s. By induction, the union of the subjacent sets of the 
blocks with minima < s is the interval s — 1. The last successive minimum of these blocks is 
< s. We insert (s) immediately after the last of these blocks. That will make ((s)) a singlet 
block. Continue until the set S is finished and one obtains a registry with range n. 

Theorem 7.2. (the stripping lemma) If TZ is a registry with support n in which the only 
singlets are singlet blocks, if one strips away these singlets and then one reintroduces them 
with the above algorithm, then the resulting registry is TZ. 

The main result of this section is the following decomposition of any permutation into an 
ordered list of very simple pieces, families, with no restrictions on the image. 

Theorem 7.3. Permutations vr o/n are in bijection with the stripped registries TZ = (Fi, . . . , Fk) 
for which, as sets, |Jj Fj C n. 

The bijection maps each singlet block of tt into the empty registry. 

It maps each nonsinglet block of tt into a registry block with the same support, with pattern 
the image of the pattern of that block. 

The images of individual blocks are concatenated in order. 

The algorithm applies to individual blocks separately, and the relative order and colors 
within the resulting families depends only on the relative order within a block. The last 
element of the block, which is its minimum, is in the last set of the corresponding block of 
the registry. 

As such, the algorithm can be extended to arbitrary injective maps which have finite 
blocks between successive minima. 

Composed with the cycle transform, the algorithm can be extended to arbitrary permu- 
tations of linearly ordered sets which have finite cycles. 

The idea of the algorithm 

We create out of a permutation, inductively by inserting the largest element, an ordered list 
of mutually disjoint families which is a registry. In it the singlet families carry no information 
and can be stripped. 

In the opposite direction of the algorithm, we map such a registry back to a permutation. 
We remove the largest element, reduce the remaining registry to a smaller registry, map that 
registry to a permutation by induction, and insert into the latter the largest element. 
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The two transforms are inverse to each other. The input of the reverse algorithm is an 
arbitrary set of mutuahy disjoint nonsinglet subsets of n, organized as famihes by specifying 
the number of ascents of each. The remaining elements of n are singlet families and can be 
added to the stripped registry before the procedure, on set places. 

In the bijection the ascents and descents (reds and blues) of the families in the registries 
end up precisely as the nonsinglet ascent respectively descent values of the permutation. 

After an inverse cycle transform the ascents and descents of families become the values 
over and respectively under the diagonal, while singlets map to values on the diagonal. This 
way the algorithm allows us to count and construct all permutations with specified elements 
under and over the diagonal. 

The decomposition of the graph of the permutation into families is well-defined and canon- 
ical. Arranging the families in a linear order, so that they form a registry, leaves a certain 
amount of freedom. There are several different algorithms which are reversible, i.e., they result 
in bijections between registries and permutations. We shall use an algorithm which, although 
more elaborate, has the advantage that it preserves in the registry much of the structure of 
the permutation itself. 

There are two distinct aspects of the algorithm, each of which is reversible. 

In the direct algorithm the highest new element of the permutation, when it does not fit 
into the family of the friend, pulls out its friend to start a new small family of two. In the 
reverse algorithm, the highest element is taken out of its small family of two, and leaves the 
friend, the other member, as an orphan in search of a family. 

The first aspect concerns the location of the newly created family and of the remaining 
old family in the registry. In the reverse algorithm, the orphan looks for and finds a family 
to enter and then places the newly enlarged family in the registry. These two processes must 
be inverse to each other. There are different possibilities for this step which lead to different 
variants of the algorithm. 

The second aspect concerns the position inside its family of the friend who is taken out to 
form a new small family. This is the content of the entry/exit lemma. The positions inside a 
family from which a friend cannot be taken out are the highest blue and the unique red. For 
the reverse algorithm, the orphaned friend enters its newfound family, and finds its relative 
position and its color in that family (see the illustration). Only the value (the height) of the 
orphan is known. The regular entry algorithm directs the orphan in its new family precisely 
on the positions other than the highest blue and the unique red. This ensures that exiting a 
family from above in the direct algorithm, and entering regularly a family from the side in 
the reverse algorithm, are inverse to each other. 

Definition Given a family F, m £ F and an element n > niaxF we say that n fits (or 
can enter) into the family F with m as its friend if inserting n immediately before m into the 
family F written as a sequence satisfies the definition of a sequence family. 

Whether n fits in F with m as its friend depends only on m, not on n. If F = (m), n 
always fits. In the entry/exit lemma, we called the positions of m for which n does not fit 
exit positions. Let us reformulate that lemma as follows. 

Given a nonsinglet family F, m € F and an element n > maxF, let F' denote the family 
F with m removed. The following are equivalent. 

n does not fit into F with m as a friend 

F is obtained by the regular entry of m into F' . 

Example 11 fits into (6,8,9,3,2) with 3 as friend, since (6,8,9,11,3,2) is a family, but 
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not with 6 as friend, since (11, 6, 8, 9, 3, 2) is not a family. The latter property is equivalent by 
the entry/exit lemma with the fact that inserting regularly 6 into its complement (8,9,3,2) 
gives (6,8,9,3,2). 



Families 



A set family 

oT 9C 



A sequence family 

9J, 



Entering a family from above 







A bicolored set. 
At least 
one red and 
one blue; 
All reds 
larger than 
all blues 



famUy 
' anchor 
(6,8,9,3,2) (smallest) 
A sequence, 
reds = ascending values 
(first or larger than previous) 
blues = descending values 
(smaller than previous) 
Reds located first. 
All reds larger than all blues. 
The order in the sequence 
determines the colors 




Entering the family The highest (new) element 
of a blue friend entered the family 

from above. As a red of its blue friend (= element 
between reds and blues located to its right.) 
(always possible) This creates an extra red. 

Exit positions from above 



Entering the family The highest (new) 



of a red friend 
from above 
at the beginning 
(when there is a 
unique red.) 



element enters the family of 
its red friend (to its right.) 
The friend changes color. 
This creates an extra blue. 



The start of a new family 

highest (new) 

ill 



highest 

(new): 



fit fit fit fit 



cofounder 

'"?'' ^', \ New family 

in old \ ' 

family9j»fo"nder 




IJJvt t 



Exiting a family 
from above: 
the highest (new) 
does not fit 
into the family 



active family 

A node of the family tree 
(its pattern is a 
primitive permutation) 



fit fit fit 

* ** * m 

' Q * *• (the crossed out 

: unique red A ; ; exit positions 

, • , . . . . \ • • are precisely 

the highest blue \ ; ; *. ^ 

ana the unque red \ : : ^he entrance 

ar^ never taken out \ • : positions for the 

V • regular entry 

highest blue J\ ; from the side 

^^^^^ shown below) 

If its friend is the highest blue or the unique red, 

the highest new element enters the family. 

Else the friend is pulled out (exits) and they form a new family. 

They form the root of a new branch of the family tree. 

In the old family they are viewed as a bud. 

The new family develops independently of the old family 

or of its sibling offshoots (other than the position of its root) 




highest blue 



Entering a family from the side 



Regular entry colors from the side 





(6,8,9,3,2) 



(6,7,8,9,3,2) 



lowest red or 

highest blue "^ yj 

highest blue 
friend here 
blue could have not 

been taken out in -^-C 
an entrance from above 



red 



red 

(not the unique red) 

blue 

(not the highest blue) 



Entering a family The new element 



from the side. 

The height 

(value) 

of the entrant 

is given. 



can only enter 
as a red. 
This creates 
an extra red. 



Entering a family 
from the side 
Possible colors 
after entry 



In the regular entry, 

the middle entrance is chosen to be red. 

This choice makes the algorithms bijective, 

since in the opposite direction, 

the entering friend will be taken out from above. 

All the positions after a regular entry 

are taken out (crossed) in the direct algorithm above. 



Definition We call the smallest element of a family or block of a permutation or block 
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of a registry the anchor. The anchor is located on the last position in a sequence family or 
permutation, and in the last position of the last family of a family block. 

Definition Call a family F in a registry 7^ in a left slide position if it is not on the first 
position and if the anchor of the family F' immediately to its left is bigger than the anchor 
of F. The name suggests that F can slide under F' when ordered by their anchors. We say 
that F has slid to its left wall if F has exchanged positions in the registry, successively, 
with families immediately to its left under which it can slide, until it is no longer in a slide 
position, i.e. is either on the first (leftmost) position, or the family to its left has a smaller 
anchor than F. 

Define similarly a right slide position for F and the slide of F to its right wall. 

Example In the registry ((9, 3, 1), (7, 4), (8), (6, 5, 2)) the family (6, 5, 2) is in a left slide po- 
sition, since 2 < 8 and slides left as ((9, 3, 1), (7, 4), (6, 5, 2), (8)) and ((9, 3, 1), (6, 5, 2), (7, 4), (8)) 
where it is at its left wall, since 1 < 2. We now start to describe the two directions of the 
algorithm. 

Permutation to Registry (the direct algorithm). 

We start with a permutation vr of n. Let vr' be the permutation obtained by removing 
from pi, written in line notation, its highest value n. We assume by induction that a registry 
TZ' was constructed for vr'. We build a registry 7^ for vr. 

Case (A). If n is last in vr then we let n form a singlet family (n) which is placed at the 
right end of the registry TZ' to form the registry TZ. Remark that Case (A) is precisely the 
case in which n forms a singlet block in vr. 

If n is not last, the element m which is positioned after n in vr is called the friend of m. 
Let F' be the family of the friend miuTl! . 

If n can be inserted into F' with m as friend, we insert n immediately before m to obtain 
a family F. We replace F' in the registry TV by F to obtain TZ. 

We distinguish two cases. 

Case (B). F' is a singlet [m) which by induction will be last in IZ' . Inserting n immediately 
before m produces a last block (n, m) of 7Z. 

Case (C). F' is not a singlet, and n can be inserted into it in front of m to form a family 
F. 

If n does not fit into the family F' with m as friend, m is pulled out of F' to obtain a 
family F, and m together with n form a new family N = (n, m). The old family F' is removed 
from TZ' and F and A^ are inserted to form the registry TZ. 

For the positions of F and N into TZ we distinguish two cases. 

Case (D). If the friend m was the anchor of its nonsinglet family F' , in 7l! replace F' by 
the pair {F,N). Slide F' (with m in it) to its left wall, then remove m from it to leave F. 
Leave the new family N in the former place of F' . This gives TZ. 

Case (E). If the friend m was not the anchor of its family F' , in TZ' replace F' by the 
pair (A, F). Slide A to its left wall and leave F in the place of F' to obtain TZ. 

This ends the direct algorithm. 

Comments 

Assume by induction that the blocks of TZ' were, in order, transforms by the above algo- 
rithm of the blocks of vr'. In Case (A) a new singlet registry block (n) was appended to TZ' 
at the end, and vr' has the same n at the end as a block. This comes from the fact that the 
last successive minimum of both vr' and TZ' is their last element n', and when appending n 
the element n' will be the next to last successive minimum, with n the last one. 
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In the other cases n, bemg largest, does not change the successive minima. In Case (B) 
the new family F stays in the place of the old one F' . In Case (D) and Case (E) the new 
families F and A^ never move to the right of the former position of F' . The family among F 
and N which contains the smallest element of F' remains in the place of F' . The other one 
(which could not contain one of the successive minima) slides to the left only past families 
with a larger minimum or anchor, so it remains within the same block. 

Thus the blocks of the registry TZ have as support the blocks of the permutation vr, in 
order. As the step of the algorithm worked only within a block and depended only on the 
relative positions and relative values of the elements in it, the algorithm transforms each 
block TT^' of TT into a block TZ^' of TZ independently of the other blocks, and the pattern of 
the resulting TZ^^' depends only on the pattern of vr''^ 

This shows that it is enough to know the algorithm on single block permutations, i.e. on 
permutations with 1 on the last position. These will be mapped to registries with one block, 
having 1 on the last position of the last family. From this, the case of a general permutation 
vr is obtained as follows. For each block vr"^ of vr, let (j)^' and ■0 denote the order preserving 
maps for which vr'-'' = ^''^ o vr''' o i/;W is a permutation, the pattern of vr'''. Let 7^'" be the 
registry produced by the algorithm from tt^'''. Let TZ^^' be the registry obtained by applying 
((j) ) to each element of every family of TZ^^'. Now concatenate "/^ = U; ^ to obtain the 
final registry. 

Since the singlet blocks of vr are mapped into the singlet blocks of TZ in order, their image 
carries no information in TZ and can be stripped off. Conversely, given a registry TZ' with no 
singlets and a range n, the complement of the range of 7?.' in n, organized as singlet families, 
can be inserted as blocks between the blocks of TZ' in order to obtain a registry 7^ in which 
the added elements are singlet blocks. 

For the reverse algorithm, registry to permutation, given a stripped registry TZ' we com- 
plete it with singlets inserted as singlet blocks as above to make its range an interval, before 
proceeding. 

Registry to Permutation (the reverse algorithm). 

We start with a registry TZ in which the singlet families form the singlet blocks of TZ. We 
construct, distinguishing cases corresponding to the ones in the direct algorithm, a registry 
TZ' not containing the highest element n from 7^. We construct by induction a permutation 
vr' from TZ' and we insert n into vr' to obtain vr. 

We distinguish four cases. 

Case (A)', n is in a singlet family (n) of TZ. By our hypothesis, (n) is a block of TZ, 
so since n is the largest element in 7^, (n) will be placed last in TZ. Remove (n) from TZ to 
obtain TZ' , map by induction TZ' to a permutation vr' and add n at the end of vr' to get the 
permutation vr, in which n will form the last block. 

In the rest of the algorithm, n is in a family F which is not a singlet. Being largest, n 
cannot be on the last place in its family. Call the friend of n the element m immediately to 
its right in F. We shall remove n, find a family for its friend m and position it in the registry 
to obtain a registry TZ' , map TZ' to a permutation vr', inductively, and then insert n in vr' in 
front of its friend m to obtain vr. 

We describe the way in which we obtain TZ' for the inductive step. 

Case (C)'. n is in a large family F of TZ. Being largest, n cannot be on the last place in 
F, written as a sequence. Remove n from F to obtain a nonsinglet family F' , and call TZ' the 
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registry obtained by replacing F by F' in 7^. 

Remark that inserting or removing the largest element n on any position other than as a 
singlet block (which case was treated as Case (A)) does not change the block structure of tt 
or TZ. 

In the remaining cases, n is in a small family N = (n, m) of TZ with its friend m. Removing 
n from N leaves the friend m as an orphan looking for a family. We distinguish two cases. 

If [m) as a singlet family placed in the position of iV in 7^ is not in a left slide position, 
then let {m) slide to its right wall. 

Case (B)'. If (m) slid right to the last position, keep it as a singlet family there to obtain 
the registry TZ' . Map Ti! inductively to a permutation vr' in which {m) will be the last block 
(and thus m is the largest in TZ' .) Insert now n before its friend m in vr' to obtain vr, where 
(n, m) will be the last block. 

Case (E)'. If {m) slid to its right wall and is not last, let F be the family to the right of 
(m). F cannot be a singlet (p), since in that case {p) would be a singlet block. A^ = {n,m), 
being to its left, would be in an earlier block, so the minimum m oi N would be smaller than 
(p) so (m) could have slid past (p) to the right, contradicting the hypothesis that (m) was at 
its right wall. So F is not a singlet family. Insert now in a regular way m into F to obtain a 
large family F' and thus obtain 7^'. 

Case (D)'. In the remaining case (m) as a singlet family placed in the position of A in 7^ 
is in a left slide position. Make (m) slide to its left wall. Let F be the family immediately to 
its right. Insert m into F in a regular way to obtain a family F'. As (m) has slid past F, the 
anchor of F is larger than m. Thus when introduced in F, m will become the new minimum, 
or anchor, of F' . Move now F' on the position occupied by A in 7?. to obtain TZ' . 

This ends the inverse algorithm. 

We now prove that the algorithm is bijective, case by case. In the direct case, n is 
the largest element of a permutation vr, which produces a registry TZ. Removing n leaves a 
permutation vr' which maps to a registry TZ' . By induction, we may assume that TZ' gives back 
vr' by the reverse algorithm. We have to show that TZ gives back vr by the reverse algorithm. 

In Case (A), the new element n is at the end of the permutation vr. In the registry TZ, n 
forms a singlet placed at the end of the registry Tl! . In reverse, this is Case (A)' in which n 
is in a singlet family (n) at the end of TZ and is removed, to obtain TZ' . That gives back vr', 
in which n is inserted at the end to obtain vr. 

In the remaining cases, n is not last in vr and the element following it is called its friend 
m. n is removed to obtain a permutation vr' which maps to a registry TZ' . We also mapped vr 
to a registry TZ. 

For the reverse algorithm, we define a friend of n in TZ. We must make sure that this 
friend is m. Then TZ is reduced to a registry not containing n. We must make sure that this 
registry is 7^'. 

By induction, TZ' is mapped back onto vr', into which n is inserted in front of its friend m, 
which gives back vr and ends this part of the proof. 

In Case (B), m is last in vr' and forms a block. Its family F' in TZ' is a singlet (m) positioned 
at the end, as a block. In TZ the last block is {n,m). Thus in the reverse algorithm, m will 
be the friend of n, and we shall be in the Case (B)' in which after removing n, the remaining 
family is the singlet (m) positioned at the end. This is TZ' , and we are done. 

In Case (C), m is in TZ' in a nonsinglet family F' , n is inserted into F' in front of m to 
form a large family F which is replacing F' in TZ' to obtain TZ. In the reverse direction n is 
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in the large family F, so we are in Case (C)'. The element following n in F is its friend, and 
is thus m. The new registry is obtained by removing n from F, and is thus Tl! , and we are 
done. 

In the remaining cases n does not fit into the family F' of m with in as a friend, m is 
pulled out leaving the family F behind, and forms a new small family N = (n,m). In the 
reverse algorithm n will be pulled away and its friend will be m which is left as an orphan. 
We have to show that the registry obtained from 11 isTZ' . 

We are in the cases (D)' and (E)', as the other cases were exhausted. In the inverse 
algorithm, n will be removed from A^ to leave its friend m an orphan. This is the most subtle 
part of the argument. 

In Case (D), the friend m was the anchor of its nonsinglet family F' . In IZ' F' was replaced 
by the pair {F, N). F' with m in it slid to its left wall, then m was removed from F' to leave 
F and joined N which stayed in the former place of F' to obtain TZ. Note that in this case, as 
m was the smallest in F' , the smallest remaining element of the remaining F will be bigger 
than m. Moreover, F will slide to its left wall past families with even bigger minima. In any 
case, the family to the left of m has a minimum larger than m, and this shows that we are in 
Case (D)', since (ttt.) can slide to the left, (m) does now slide past families with larger minima 
to its left wall, which is precisely where it left the rest F of its former family F' . (m) makes 
one last left slide past its former family F which has a minimum larger than m, before ending 
at its left wall. F is now immediately at m's right, where in Case (D)' m finds it. Now m 
enters i^ in a regular way. The entry/exit lemma guarantees that m will end up precisely in 
the place from which it was taken out by n in Case (D), in this case the anchor, or smallest 
position, and thus the family F' is reconstituted. Then the family F' is moved on the position 
of A^ = {n,m) in R which is precisely the place from which it was taken away in Case (D). 
So with F' , as it was before, in its old place we recovered 7^' and we are done. 

In this case the idea was that the family head m, in a position of responsibility, could not 
leave home. So it sent its family on vacation South, to the left (block anchors get smaller, 
and thus move Southward to the left). To make sure that its family doesn't get lost, m 
accompanied them all the way to ttt-'s left wall and left them there, where they would be easy 
to find. Back home, m remembered the situation by the fact that there was slide room to its 
left. Home alone, m started a new family with n, but that didn't go very far, since n left. So 
m went back, sliding all the way to its left wall, to find its former family there. Its former 
family gladly gave m the honor position which it used to have. Then m brought them all 
back to their old homestead, and all was happy like before (which is what bijections are all 
about.) 

In Case (E), the friend m was not the anchor of its nonsinglet family F' . In this case 
N = {n,m) slid to its left wall. Thus after n is removed, (m) cannot slide further left, which 
shows that we are in Case (E)'. So (m) slides right. It moves to its right past all the families 
under which (n, m) slid left, but cannot slide right past its former family F, since m was not 
the anchor, so F has a smaller minimum than m. Thus m has found its former family F at 
the right of its right wall, and enters it regularly. By the entry/exit lemma, m recovers its 
former position, to form F' in the place of F, which is precisely the place where F' used to 
be in 7^', and we are done. 

In this case the idea was that a family member m, not being the family head, took some 
time off and went on vacation to its left wall, m remembered that by the fact that there was 
no more room to go to further to its left. While there, m started a family with n, but as 
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vacation flings go, n left. So m went back sliding to its right, to be stopped only by its former 
family. It reentered it in its old place and, once again, all was the way it used to be. A happy 
ending for bijective proofs. 

This ends the bijectivity proof in one direction. 

For the other direction, we start with a registry TZ. We remove from it the largest element 
n and change it to a registry TZ' . By induction 7^' is mapped onto a permutation vr' which 
by the direct algorithm gives back TZ' . When n is inserted into vr' to give vr, and vr is mapped 
onto a registry, we must show that that registry is TZ. 

In Case (A)', n is in a singlet family (n), the last family of 7^. 7^' is obtained by removing 
(n), is mapped into the permutation vr'. n is inserted at the end of vr', so we are in Case (A). 
The new registry is obtained by appending (n) at the end of TZ' , and that is precisely TZ. 

In the remaining cases, n is in a family F (written as a sequence) which is not a singlet, 
and has a friend m to its right in i^. n is removed from F to obtain a new registry TZ' which 
is mapped into the permutation vr'. n is introduced in vr' immediately before its friend m to 
obtain vr. By induction, vr' is mapped back into 7?.'. We need to show that vr is mapped into 
7^. 

In Case (C)', the family F of n is large, with > 3 elements, n is removed from it to get 
a family F' with > 2 elements, which replaces F and thus one obtains 7^'. So F' is not a 
singlet. Inserting n back in front of its friend m produces F, which is a family. So n can be 
inserted, and we are in Case (C). We replace F' with F and we obtain the registry 7^ from 
which we started. 

In the remaining cases, n forms, with its friend m, a small family N = (n, m) of TZ. 
Removing n from A^ leaves the friend m an orphan. 

If (m) as a singlet family placed in the position of A^ in 7^ was not in a left slide position, 
we let (m) slide to its right wall. 

Case (B)'. If (m) was on the last position, we kept it as a singlet family there to obtain 
the registry 7^'. We mapped TZ' inductively to a permutation vr' in which (m) will be the last 
block (and thus m is the largest in TZ' .) We inserted n before its friend tti in vr', to obtain vr 
where {n,m) will be the last block. This is precisely the Case (B). vr' is mapped to TZ' , and 
the registry obtained by replacing (ttt.) by (re, rre) in TZ' is precisely the 7^ we started from, so 
vr is mapped into TZ. 

Again, the remaining cases are more delicate. We started from TZ in which n was in a 
small family A^ = (re, m). We removed re and were left with (rre). 

Case (E)'. (?re) could not slide left, so (m) slid to its right wall and was stopped by a 
family F from sliding further right. F was a nonsinglet family into which rre was inserted in 
a regular way, to obtain a large family F' . With F' in the place of F, we called the registry 
TZ'. 

Note that to get back TZ out of TZ' , we pull m out of F', leave the remaining F in place of 
F' , start a small family (re,rre)and slide it to its left wall. 

TZ' was inductively mapped to a permutation vr' which was then sent inductively back to 
TZ' . We have to show that the permutation vr obtained by inserting n in front of its friend rre 
is mapped by the direct algorithm back into TZ. 

Since rre was inserted in a regular way into F to get F' , by the entry /exit lemma rre is in 
an exit position in F' , that is, when re is placed in front of m in F' , re will not fit and will pull 
its friend rre out of its family to leave back F. In TZ' , (rre) could not slide past F, so rre was 
larger than the minimum of F. So m was not the anchor of F' , which places us in Case (E). 
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In Case (E), F' had m pulled out by n. So after pulling m out of F' we get back F. 
In Case (E) A^ = {n,m) was slid to its left wall (where it cannot slide further left), and F 
replaced F' to get a new registry, so sliding back m to its right wall it finds F and enters 
it regularly, by the entry/exit lemma on its former position in F' . These operations are 
precisely the reciprocal for Case (E)', 

From our original TZ, in Case (E)' we obtained 7^' from 7^ by removing n from its small 
family {n,m), by having the orphaned m, which could not slide left, slide right to its right 
wall, in front of a family F into which m was inserted regularly to get a family F' . Thus 
reciprocally, to get TZ from 7^', we remove m from F' to leave F, form the new family (n, m), 
and slide it to the left wall where it cannot slide further left, and from where, when (m) is 
sliding right, it will be stopped by F. These are precisely the transformations in Case (E), so 
those lead back to our original registry TZ. 

Case (D)'. Here {m) as a singlet family placed in the position of A^ in 7^ is in a left slide 
position. Make (m) slide to its left wall. Let F be the family immediately to its right. Insert 
m into F in a regular way to obtain a family F' . As (m) has slid past F, the anchor of F is 
larger than m. Thus when introduced in F, m will become the new minimum, or anchor, of 
F' . Move now F' on the position occupied by A^ in 7^ to obtain TZ' . 

The reverse of these operations is the following, giving back TZ out of TZ' . Take the family 
F' and slide it to its left wall. Remove m from it to remain with F. Slide m to the former 
position of F' and insert n to make the family (n, m) there. This gives back 7?.. 

As m was the anchor of F' and n did not fit before its friend m into F' (since m had 
entered regularly into F to get F') we are in Case (D). The operation in that case for getting 
TZ out of 7^' are precisely the ones we described earlier as reversing Case (D)' to get the 
original TZ out of TZ' . So Case (D) gives back our original TZ and we are finished. 

This ends the proof of the bijectivity of the algorithm. 

The algorithm put arbitrary permutations of n in bijection with registries in which singlets 
were blocks. By the stripping lemma, the latter are in bijection with arbitrary stripped 
registries, and the main theorem is proved. 

The pattern of a nonsinglet permutation block is a permutation of p which ends in 1. By 
the inverse cycle transform, this corresponds to a nonsinglet single cycle permutation. The 
algorithm, as we described it, maps single blocks into single block registries, i.e. registries 
having 1 in the last family. The number of such registries, for a given number total number 
k of reds, or ascent values and I of blues, or descent values, is the shifted multinomial A'^^ ^ , 
since all families except for the one containing 1 can be positioned arbitrarily. The shifted 
multinomials A^^ ^ are the Eulerian numbers E^^i- 

In general the algorithm as given works on block patterns. That is, in order to map an 
arbitrary permutation vr, decompose vr into blocks vr'". Bring each block, of size rii, to its 
pattern, which is a permutation of nj which ends in 1. Lift back the resulting registry TZ^^' to 
the original values taken by vr^", and then concatenate the results in order. 

In what follows we write permutations in line notation, though written enclosed in paren- 
theses. We define inductively a primitive permutation 

Definition A permutation vr of n is called primitive and has active family ^(vr) if 

n = 2 and vr = ^(vr) = (2, 1) 

or n > 2 and denoting by vr' obtained from ir by removing n, 

the permutation vr' is primitive, and 
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n is followed in vr by tti E -^{t^') its friend, and either 

Case (A) ^(vr') with n inserted (in the order from vr) is a family, which is then A{Tr), or 

Case (B) ^(vr') with n inserted is not a family, and then ^(vr) is ^(vr') with m removed. 

In the latter case, (n, m) is a small family called a bud of vr. The buds of vr are the buds 
obtained at all the steps, when building vr out of (2, 1). The buds are mutually disjoint, and 
their complement in vr is the active family ^(vr). 

The pairing pattern 'P(vr) of the primitive permutation vr is a sequence of integers of 
length n, defined inductively as follows. 

If vr = (2, 1) then 'P(vr) = (1, 2), and if n > 2, let b denote the number of buds, and / the 
length of .A(7r), so 2b + f = n 

in Case (A) the pattern 'P(vr) is obtained by appending 6 + / to 'P(vr'). in Case (B) 
the pattern 'P(vr) is obtained by appending to 'P(vr') the negative number —V{'K'){m) (where 
V{TT'){m) is the m-ih. entry of V{tt').) 

In Case (A), the last entry in V{'k) is positive, not paired, on the n-th place with n part 
of the active family ^(vr). In Case (B), if —p is the last entry in ^(vr), the pair {p,—p) in 
V{tt) is on positions {m,n) and thus we have. 

Define a pairing sequence V = {pi,p2, ■ ■ ■ ,Pn) as follows. 

The first 2 entries of V are (1,2). 

The positive numbers appear in natural increasing order {1,2,...} in P, 

Each negative number is paired with a previous positive number, with different negative 
numbers having different pairs. 

The sum of the signs of the entries on positions 3, . . . , A; of "P is nonnegative, for all k. 

Let sgn = sgn{V) denote the sequence of signs of V. For each descent sgn[i) = —1 there 
are Ylij<i ^d'^^iJ) choices for Vii). 

We may think thus of {sgn[2>),sgn{4), . . .) as a Dyck path having nonnegative partial 
sums, with a decoration of any descent by a number between 1 and 2 + the path height 
before the descent. 

Remark now that the pairing sequence V{tt) of a primitive permutation vr satisfies the 
above definition of a pairing sequence. 

The sum of the signs of the numbers is the size / of the active family ^(vr), and we always 
have f >2, as f = k + l with k and I the number of reds and blues (ascent values and descent 
values) of w4(vr), and k,l > 1 

The positions of the unpaired positive numbers in 'P(vr) are the numbers in the active 
family ^(vr). 

The number of negative numbers in 'P(vr) is the number of buds b. 

The positions of each pair (i, —i) in P(vr) form a bud of vr. 

Our main result is the following. 

Theorem 7.4. Pairs {V, {k, I)) between a pairing sequence V and numbers k,l > 1 with k-\-l 
equal to the sum of signs of V are in natural bijection with primitive permutations. 

The bijection matches {V, {k,l)) with a primitive permutation vr having pattern V{tt) = V 
and with an active family .A(vr) having k reds and I blues. 

It is implemented by an algorithm which constructs the primitive permutations directly, 
without going through all permutations. 

The singlet families in the registry thus carry no information other than the subset which 
consists of their union. The algorithm establishes a bijection between permutations and 
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registries consisting of families with > 2 elements, i.e., mutually disjoint subsets with > 2 
elements, each with a specified, nonzero number of ascents and descents. Singlet families are 
inserted into such a registry in a canonical way, with each singlet starting with the smallest 
one inserted immediately to the right of the set of families in the registry with minima smaller 
than it. 

Remark that the algorithm described above maps each block of the permutation into a 
block of the registry, independently of all the rest. Thus the algorithm applies to permutations 
of the integers Z which have finite support, or by the inverse cycle transform, to permutations 
with a finite number of cycles. More generally, the algorithm applies to permutations of Z 
with finite cycles, by transforming each cycle individually into a block of the registry. 

An example. In the sequel, we mark ascents in red and descents in blue for clarity. We 
start with the permutation 6251473, with 567 over diagonal in red, 2 on the diagonal in black 
and 134 under the diagonal in blue. Its cycle transform is tt = 6735412, with 567 now ascent 
values in red, 2 a singlet block in black and 134 descent values in blue. We work with the 
cutoff vr^^ of vr to values < k, for k = 1,2,..., and think, following Euler, of vr^^ as obtained 
from vr^'^"^) by inserting the largest element k. We construct the families and the registries 
(lists of families) TV- ' corresponding to each vr' •* . 

Start with the singlet vr'^^ = 1 so Tiy' = ((1)). The next element is 2, with permutation 
7r(2) = 12. As 2 is last, it is inserted in the registry as the singlet (2) at the end, so TZ^"^' 
= ((1),(2)). The next cutoff is vr'^^ = 312. The friend of the largest element 3 is its right 
neighbor 1, with family (1). 3 joins it and gives the family (31), remaining in the place of 
(1), in the first place in the registry, which now becomes TO-'^' = ((31), (2)). Then -k^^' = 3412. 
The friend of the largest element 4 is its right neighbor 1, with family (31). With 4 inserted, 
that family becomes (341), which is a family. So TZ^^' = ((341), (2)). 
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The next step is vr^^^ = 35412. The largest element 5 is followed by the friend 4. When 
inserted into the family (341) of its friend, we obtain (3541). This is not a family, as the 
descent 4 is bigger than the ascent 3. So 5 pulls its friend 4 out and they start a new family 
54. As the element 4, taken away, is not the anchor of or minimum of its family, that family 
remains in place and the new family (5, 4) is sliding to its left wall, i.e., until it can move no 
further under elements with an anchor bigger than 4. In this case, the left wall is the left end 
of the registry, and (5, 4) ends up in the first place in the registry. Note that the idea behind 
this choice of position is the following. In principle, in the reverse algorithm, 4 will look for its 
old family, which could be to its left or to its right. The fact that (5,4) cannot slide further 
left is a giveaway for the fact that it was not an anchor of its family, and that it must slide 
back right, stopped by its old family, which it will then rejoin. 
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Thus, the new family (5, 4) is placed first in the registry, to give TZ^^> = ((54), (31), (2)). In 
the next step, vr^*') = 635412, 6 enters the family of its friend 3, with 7^(^) = ((54), (631), (2)). 
In the last step vr = tt^^' = 6735412, 7 enters the family of its friend 3, with TZ = 7ly' = 
((54), (6731), (2)). 
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Each step is reversible, so the permutation can be reconstructed from its registry 

For the reverse algorithm, we construct out of each registry a smaller registry, by removing 
its largest element, and which we map to a permutation by induction. Then we insert the 
largest element into that permutation. 

We check the steps as follows. The registry TV-"^' = ((1)) gives the permutation vr'-"'^^ = 1. 
^ = ((1), (2)) gives the previous registry Tiy' = ((1)), and then vr^^^ = 12 by inserting the 
singlet 2 at the end of vr^^^ The registry W-^' = ((31), (2)) leaves, after removing the largest 
element 3, its orphaned friend 1 which succeeded it in the family (31). The orphan is left in 
place, so TZ^'^' = ((1), (2)). By our previous step, inductively, IZ^"^' corresponds to vr'^^ = 12. 
The largest 3 is placed in the permutation vr'^' = 12 immediately preceding its friend 1, to 
obtain 7r(3) = 312. Now 7^('^) = ((341), (2)) has largest element 4 followed by its friend 1. 
With 4 removed from its (large) family (341) we get 7^^"^^ = ((31), (2)) which by induction, 
as we saw above, gives the permutation tt^'^' = 312. In tt^^', 4 is inserted before its friend 
1 to get vr^"^) = 3412. Take now 7^(^) = ((54), (31), (2)). Its largest element 5 has friend 4. 
Removing 5 leaves 4 orphaned. 4 cannot slide further left, so it will slide right to find its 
family as (31), which it enters, to give (431), with 3 changing color from red to blue. The 
registry is now 'RS'^' = ((341), (2)). By induction we showed that it corresponds to vr'-^) = 3412, 
and we insert into it 5 before its friend 4, to get vr^^) = 35412. In 7^('^) = ((54), (631), (2)) 
the largest 6 has friend 3. With 6 removed we get TZ^^' = ((54), (31), (2)) which we showed 
corresponded to vr*^^^ = 35412. Inserting 6 before its friend 3 in vr'^^ gives vr'^^ = 635412. 
Finally TZ = 7^*^^^ = ((54), (6731), (2)) has its largest element 7 in a large family (6731) 
having as friend its successor 3. With 7 removed, the registry is TZ^^' = ((54), (631), (2)) 
corresponding by induction, as we saw, to vr'^-* = 635412. By inserting 7 in front of its friend 
3, we obtain the final result of the inverse algorithm, the permutation tt = vr^^^ = 6735412, 
which corresponds to the registry 7Z = ((54), (6731), (2)). 

The inverse algorithm starts in fact backwards, with the registry TZ = TZ^'^' and proceeds to 
define the registries TZ^^' , . . . , TV-^' inductively as above, and then constructs tt'^-* , • • • , t^^"^' = vr 
as described. 
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At each step k the ascent/descent colors in vr^*^^ match the colors within the families in 

The inverse algorithm thus reverses the map from vr to 7^ from the direct algorithm. 

A family registry has singlets, which are elements of the permutation which form a block 
of 1 between successive minima, as singlet blocks of the registry, again as blocks of 1 be- 
tween successive minima of the registry, as such their position is entirely determined by the 
nontrivial families in the registry. If we want to have a bijection between permutations and 
unrestricted registries, we have to strip these singlets away from these unrestricted registries. 
The reverse algorithm requires that singlets be present in the correct places, but as their 
positions are determined by the other families, the missing singlets can be inserted in place 
at the beginning of the algorithm. A stripping and reinsertion algorithm for singlets is part 
of the implementation of the main theorem. We formalize it in this paper as the stripping 
lemma. 

Essentially the algorithm decomposes a permutation vr of n into an arbitrary 
ordered list {{Fi,ai), . . . ,{Fk,ak)) ^vhere Fi,...,Fk C n are mutually disjoint and 
1 < flj < \Fi\ for all i. The complement of IJ^ Fi consists of singlets. The largest Oj elements 
in each Fi are its ascents, the rest its descents. The algorithm maps these ascents and descents 
of families into the ascent-values ttj > ttj-i respectively descent-values ttj < 7rj_i of vr, with 
singlets becoming the singlet blocks. In the inverse cycle transform o" of vr these become 
elements Uj with Uj ^ i and Ci = i respectively. The registries can be easily counted, resulting 
in the multinomial type formulae given at the beginning. 

8 The parenthesized form and the tree form of a permutation 

The algorithm presented in the previous section reveals information on the internal structure 
of a permutation, decomposing it into a tree. Every node of the tree is called a primitive 
permutation. 

A permutation decomposed into families, as seen in the illustrations, has a two dimensional 
structure, with families positioned vertically, each above an opening in a family below, in a 
tree structure. This structure can be captured in a very natural manner, by grouping the 
permutation into nested parentheses. 

We construct the parenthesized form of a permutation vr as follows. The outermost paren- 
theses enclose the blocks of the permutation, which by inverse cycle transform correspond to 
cycles. Then, every time a new family is started in the algorithm, that family starts a new 
tree branch, and we enclose that new family in parentheses. 

The result of this algorithm is the parenthesized form of the permutation vr. 

Remark that, while the main bijection depends on a choice of algorithm for ordering new 
families in a linear order in the registry, the parenthesized form of the permutation, is, by 
contrast, canonical, as are the tree and the primitive permutation structures which are derived 
from it. 

In our previous example, the parenthesized form of vr = 6735412 is constructed stepwise 
as follows. We start with the singlet vr^i) = (1), then tt^^) = (1)(2), then tt^^) = (31) (2), 
TT^"^) = (341) (2). As 5 starts a new family with its friend 4, we have tt^^) = (3(54)1) (2). 
Afterwards vr^^) = (63(54)1) (2) and vr = vr(^) = (673(54)1) (2). 

The families in the registry TZ = ((54), (6731), (2)) are now visible in the body of the 
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permutation vr = (673(54)1)(2), as the content of each parenthesis from which the inner 
parentheses and their content have been removed. 

The outermost parentheses are the blocks of the permutation vr, corresponding to the 
cycles of the inverse cycle transform of vr. 

The figure below describes the graph of a permutation on the left and its cycle transform 
in the center. Both are decomposed into families with elements over diagonal, respectively 
ascents, in red and elements under diagonal, respectively descents, in blue. Families are 
ordered in the registry on the right. The family tree associated to the permutation is placed 
underneath. The permutation in the center is written in parenthesized form, reflecting the 
tree structure underneath. The tree structure is also shown on the graph of that permutation, 
in the background. Roots and buds of each node are marked in green with branches marked 
in brown. In the graph, the active family of each node is connected by a light purple line to 
that node. 




7 1 5 10 2 4 8 13 6 12 3 11 14 9 (7 8 (13 14 9) 6 4 (10(12 11) 3)(5 2) 1) 




9 The multiparenthesized form of a permutation 

Given a permutation vr, the main algorithm constructs a registry of families, a list of nonsinglet 
families, arranged in linear order, which encode the permutation vr. We can consider now 
the permutation vr^^^ which encodes the order of the families in the registry, relative to the 
standard order in which the minima of the families are increasing. In turn, vr'^^ can be 
subjected to the same algorithm, will produce a registry, and a permutation vr'^^ encoding 
the order of the families of vr'^-* in the registry. We can continue this way until we reach a 
family which consists only of singlets. 

It is clear that the previous sequence of registries encodes the original permutation vr. 
We do not need, however, the order of families in any of the registries, since that is encoded 
inductively in the next permutations. The only information we need to keep track of is the 
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number of reds and blues in each family, in all iterations. We obtain this way the following 
data, which encodes, bijectively, any permutation n of the set n = {1, . . . , n}. 

Partition the set n into nonempty subsets. For each nonsinglet family, specify the number 
of reds and blues (e.g., by coloring a number of its elements starting from the top as red, with 
the rest blue). 

Partition the set of nonsinglet sets from the previous step. Specify for each nonsinglet at 
this step a number of reds and blues (e.g. by coloring some of its members red, starting with 
the biggest ones, as arranged in increasing order of the smallest number in each member). We 
mark this second partition, by enclosing the previous partition in parentheses, and we mark 
the colors by coloring these parentheses. 

At each subsequent step, partition the nonsinglets of the previous partition, enclose each 
nonsinglet part in a parenthesis, and color the immediate inner parenthesis within that red 
or blue, as before. Repeat this until all elements of the last partition are singlets. This will 
encode the original permutation tt. 

The idea is that if we put the elements of a set into bags with > 2 elements each, leaving 
out some singlets, and we mark on each bag a number of reds, then we put bags into bigger 
bags, each with > 2 bags inside, leaving aside some singlets, and mark on each bigger bag 
a number of reds, and we continue this way until there are only singlets left, the data in 
the successive partitions, together with the number of reds on each bag, encode precisely a 
permutation of the original set. 

Example. For the identity permutation, the set consists of singlets and no bags are 
needed. If the permutation is the reverse, {n,n — !,...,!), then we need one bag, with one 
red and the others blue. 

In the picture of the large permutation attached in the appendix, the idea of the multi- 
parenthesized form is to start with the permutation on the upper middle board, and continue 
with the permutation on the lower right board, describing the order of the families, and then 
proceed from the latter inductively, until the lower right board becomes trivial, i.e., there are 
< 1 nonsinglet families. Thus, the main idea of the multiparenthesized form of a permutation 
is that by marking the number of reds and blues of families in successive registries, we no 
longer need to specify the relative order in any of the registries. In other words, with our 
algorithm, repeated bicolored partitions encode the relative order of families in registries, and 
thus ultimately the original permutation itself. 

11198175 12 1373918101516126 14 204 
{{11198{175))12 13 73{9{18 10)15 161)) (2) (6 14 204) 
{ { { 11 19 8) (17 5 ) ) ( ( 12 13 7 3 ) (18 10 ) ( 9 15 16 1) ) ) (2 ) (6 14 20 4) 
A permutation with its ascent and descent values colored 
The parenthesized form of the permutation 
The multiparenthesized form of the permutation 

We have implemented in Mathematica both directions of the algorithm, as well as supple- 
mentary aspects of the structure, fully interactively. One can input permutations or registries 
by clicking and dragging any point in the display. 

The tree form of a permutation. The parentheses give a permutation a quite complex 
tree structure, visible in its graph. We now construct from each parenthesis a node of a tree. 
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We shall retain at a node enough information for reconstituting the tree, and thus the whole 
permutation, and call its pattern a primitive permutation. 

A primitive permutation at a node has a branch root in a previous node, it contains a 
family, as defined before, called the active family of the node, as well as pairs of elements 
forming small families called buds, which are the branch roots of its offshoot branches. The 
idea is that while a permutation decomposes into families as seen in the algorithm, the branch 
root and buds indicate precisely where different families are inserted into previous families. 

In the example in the picture above a node is (7 8 (13 9) 6 4 (10 3) (5 2) 1), where only 
the smallest two of each inner parenthesis were retained as a bud which started that inner 
parenthesis. Its pattern is itself, so this is a primitive permutation. Its active family is 7 8 
6 4 1 and its root is 2 1. Note that in the tree underneath the permutation, in the picture 
above, the part of the permutation at each node has been reduced to its pattern, which is a 
primitive permutation. For instance the node (10 (12 11) 3) has pattern (2 (4 3) 1). 

The tree can be reconstituted from these primitive permutations by specifying for each 
branch root the bud on the previous root with which it will be identified (unless it is a root 
of a whole tree.) In addition one has to specify the shuffling of various branches, the way 
in which the numbers on different branches are made different while preserving the order on 
each branch. 

Note that the parenthesized form of a permutation is natural and canonical. There are 
various algorithms for ordering the families in the registry, which are based on the history of 
the families, as in the previous example. The "visual" orders, based on the position of the 
smallest or largest ascent or descent of each family in n do not lead to bijective algorithms. 

The structure described, as well as the definition of primitive permutations below, are 
motivated by the following property. The pattern structures of different nodes (an 
active family, buds and a root for each) are independent of each other. In other words, 
once a new family is started, it has the potential to develop like any self standing permutation. 
The only relative interaction with the rest is that, vr being a permutation, each entry can be 
used only once, i.e. different branched need to be shuffled. This shuffling, which is a pattern 
on the union of the branches which reduces to the pattern of each node, is arbitrary. 

10 Primitive permutations 

The idea of the primitive permutation is the following. In the evolution of a permutation 
which we described before, when a new element does not fit into a family, it picks its friend 
immediately to its right and they start a new family. Keeping track of all these new families 
is complicated, and as it turns out, once their new family is started they do not interfere 
anymore with their old family or with their siblings. Therefore, just like in a large modern 
family, we shall keep on the wall a wedding picture of each member who left, but we shall no 
longer follow in detail what happened to them. Such a permutation, in which a departed pair 
of immediate neighbors, called a bud (since they started a new branch of the family tree) is 
no longer developed, will be called a primitive permutation. It consists of the active family 
members still in the household, together with the positions of the buds. It forms a node of 
the family tree, with the buds indicating the precise insertion points of the branches above. 

The grand assembly of all these nodes into the family tree ~ in fact a family forest, as 
blocks of the permutation develop into separate trees - is axiomatized and described by a 
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decimal code, similar to the one which organizes a whole library. 

The delicate algorithm which describes the construction of a primitive permutation, given 
only the pairing between the values in each bud and the number of colors in the active family, 
is described in what follows. The assembly algorithm for the decimal code of a permutation 
is described in the chapter which follows. 

The data from which we construct a primitive permutation is of the following type. Con- 
struct a primitive permutation in which the fifth highest element was taken away by the 
eighth highest, and from the remaining six members the highest four are ascent images and 
the lowest are two descent images. We shall show that there is precisely one primitive family 
satisfying these conditions, and we shall construct algorithmically its whole history. 

In the decimal notation which encodes the whole permutation, this family will have a 
decimal head like 3.1.7, and then each member will be assigned two digits. In our case, those 
will be l.b, 2.b., 3.r, 4.r, 5.1, as the fifth highest was taken away, 6.r, 7.r, 5.2, which is 
on the position of the highest element which took 5 away. When 5.1 and 5.2 start their 
new branch of the family tree, the family will append to the label 3.1.7 of the family which 
they left the suffix 5, which will give the label 3.1.7.5 for their new branch. In it, they will 
be the elements 3.1.7.5.1.x, and 3. 1.7. 5. 2. y. For instance, 3. 1.7. 5. l.b will indicate that the 
founder of the new family stayed in, as a blue descent value, while 3.1.7.5.2.1 will show that 
the cofounder of the new family moved on, to start a new family of its own. 
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Primitive permutations. Recall that we named a permutation vr a family if it is a 
singlet or has the pattern d + 1, d + 2, . . . , d + a, d, (i — 1, . . . , 1 with a, d > 1 (e.g. 456321 for 
a=d=3.) Any nonempty subset of a family is a family. The formal definition of a primitive 
permutation is the following. 

We call a permutation vr of n primitive if it has > 2 elements, there are mutually disjoint 
adjacent pairs (i,i + l)ie7 in n, called buds with 7r(i) > 7r(i + 1) such that (a) on the 
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complement of the buds vr is a (nonempty) family (b) with one or more of the buds left in it, 
vr is not a family and (c) the properties (a) and (b) hold for any cutoff of vr by c = 1, . . . , n, 
defined as the pattern of the restriction of vr to points i with 7r(i) < c. The complement of 
the buds will be called the active family of the node. The smallest and the second smallest 
elements of the primitive permutation (which includes the buds) are called its root. These are 
the two elements which started the node, and will also be called the founder and cofounder. 

Primitive permutations are easy to construct inductively. Given a primitive permutation 
of n, insert n + 1 immediately before any element m of the active family. The result is always 
a primitive permutation. Either the active family together with n + 1 forms a family or else, 
if they don't, n + 1 and m form an extra bud. All primitive permutations are constructed this 
way from the singlet (1). 

The graph generating the primitive permutations is shown below. Any primitive 
permutation is uniquely built by a walk on this graph. 

Each downward edge produces family growth, by an ascent (red) or descent (blue). Each 
upward edge produces a new (parenthesized) bud and decreases the family by an ascent (red) 
or descent (blue), in a number of ways marked by the arrow number. The pattern of the 
active family at a node, with the number of ascents in red and descents in blue, is shown in 
examples. The primitive permutation may contain previous buds. In the examples, as the 
new element 6 enters before one of the 5 positions in the active family, the family may grow 
or bud. 

11 The algorithmic construction of primitive permutations 

Definition. A family pairing data is a pair (inds, URcds)- inds is a sequence with elements 
the numbers 1, . . . ,n and a subset { — !,..., —n}. The positive numbers appear in the se- 
quence in their natural order, and each negative number —k appears after the positive k. Let 
(si, S2, . . . , ) denote the signs of elements in the sequence, with Si € {±1}- For any k > 2, we 
require that the partial sum si + • • • + Sfc > 2. if s = ^ Sj, then 1 < npj,cds < s — 1. 

Remark that the sign condition is equivalent to the fact that si = S2 = +1 and (53, S4, . . .) 
form the increments of a Dyck path. Equivalently, the entire sequence of signs {si,S2,- ■ ■) 
forms the increments of a Dyck path, which reaches height 2 in 2 steps, and remains at or 
above height 2 afterwards. The number of choices in the position set corresponding to a 
descent, Sj = — 1, is precisely the height si + S2 + • • • + Sj_i of the Dyck path, before Sj. 

Thus the family pairing data is equivalent to a Dyck path which reaches height 2 in 2 steps, 
and remains at or above height 2 afterwards, in which each descent is marked with a number 
1,2, ... ,h, where h is the height of the path before the descent, together with a number 
n-Reds between 1 and the final height minus 1. An example is the Dyck path (si, . . . ,sj) = 
(+1, +1, +1, -1, +1, +1, -1), decorated as (+1, +1, +1, -1(2), +1, +1, -l(4))(nReds = 2). We 
call the preceding data a decorated Dyck path. The paired values corresponding to this data 
is the sequence inds= (1, 2, 3, —2, 4, 5, —5), where -2 is the second choice among the preceding 
(1, 2, 3), and -5 is the fourth choice among the remaining (1, 3, 4, 5), as 2 was paired with -2. 
The height of the Dyck path before the two negatives was 3, respectively 4. The final height 
of the Dyck path is h=3; we have n^eds = 2 < h, and we let reeiues = h — n^eds- 

The algorithm which follows will establish a bijection between paired values - or equiv- 
alently decorated Dyck paths - as above, having m descents, ri^eds j ^Biucs (and thus a total 
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length n = n^eds + "Blues + 27Ti), and primitive permutations with active family n^eds reds 
"Blues blues, and m buds. 

With the algorithm, the data in the previous example will correspond to the primitive 
permutation ((4 2) 3 5 (7 6) 1), with buds (4 2) and (7 6), and active family (3 5 1), having 
ascent values (reds) 3,5 and descent value (blue) 1. 

11.1 The pairing sequence — primitive permutation algorithm 

We start with paired values data as described above. The only such data of length 2 is (1,2) 
with n-Reds = "Blues = 1- The corresponding primitive permutation is vr = (2 1) with no buds. 
Assume now that paired values inds of length n and riReds are given. We let inds' denote the 
sequence obtained by removing the last entry of inds. We compute below a modified "p^^^jg 
to be paired with inds'. We let vr' denote the permutation constructed inductively by the 
algorithm from the data (inds', "j^g^s)' ^^'^ ^^ describe the construction of vr from vr'. 

Before we begin the algorithm, let us make the following observation. If the last sign 
of inds is -1, then we compute the number nchoice of its chosen sequence as follows. In the 
unpaired positive values preceding -k, its pair +k appears on position ^choice- Equivalently, 
"choice was the decoration of the last entry -1 of the Dyck path (• • • , -1 (nchoice))- In the 
example preceding the proof, (+1, +1, +1, —1(2), +1, +1, —1(4)), the last nchoice = 4. 

We shall construct n^^^^ for the previous vr' by adding a correction to n^eds- In the second 
part of the algorithm, we shall show that the permutation vr which we construct has precisely 
"Reds active reds, and thus the correction was correct. 

The previous "'j^cds ^^ computed by adding to n^eds the following correction, c = n'j^^^^^ — 

"Reds- 

Case (A) If the last sign in inds is +1 and n^cds = 1) then c = 0. 

Case (B) If the last sign in inds is +1 and n^eds > li then c = — 1. 

Case (C) If the last sign in inds is —1, and if the last choice nchoico as described above, 
is < nBiucsi the correction c = 0. 

Case (D) If the last sign in inds is —1, and nchoice > neiues) then c = +1. 

This completes the first part of the algorithm. We now have the data inds', obtained by 
removing the last entry of inds, and n^^^^^, obtained by the correction in the cases above. 

With these we construct inductively a primitive permutation vr'. We shall describe below 
the way in which we obtain the desired permutation vr from vr', and check that in each case 
the correction between n^eds of vr an n^eds of vr' is the one described above. 

Recall that n is the length of the sequence inds, and also the length of and the highest 
number in the permutation vr which we are constructing. 

Case (A) If the last sign in inds is +1 and n^eds = Ij then we insert n immediately before 
the first element in the active family of vr' (obtained from vr' by removing its parenthesized 
buds). In vr, n will be red and the previous unique red in the active family of vr' will turn blue. 
Thus, the active family of vr will have a unique red, as we assumed in the above condition 
"Reds = 1- Thus, nRcds = "Reds = 1' ^nd c = n'^^^^ - n^eds = 0, as assumed. 

Case (B) If the last sign in inds is +1 and n^eds > Ij then we insert n immediately 
before the highest blue in the active family of vr'. Then the active part of vr (obtained by 
removing its parenthesized buds) is a family, with n as red. Thus, n^eds = "Reds + 1 > 1) and 
c = "Reds ~ "Reds = -1, as assumed. 
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Case (C) If the last sign in inds is —1, with the choice number richoice a-s before, and if 
f^choice < '^Biues, let 6 denote the entry number nchoicc on the increasing list of blues in the 
active family of vr'. We replace b by the bud (n, b) to obtain the permutation vr. The number 
of active blues decreased by one, the number of active reds stayed the same, so ?iRcds 



n 



Reds' 



and c 



n 



Reds 



'^Reds = 0, as assumed. 



Case (D) If the last sign in inds is -1, and nchoicc > "-Blues, then let n[-^^^^^ = nchoice - 
"Blues > 0, and let r denote the entry number nchoice o^ ^^^ increasing list of reds in the active 
family of vr'. We replace r by the bud {n,r) to obtain the permutation vr. The number of 
active reds decreased by one, the number of active blues stayed the same, so nRods = "-Reds ~ -*■' 
and c = n'pj^eds — n^cds = +1, as assumed. 

This ends the algorithm. 

Remark that when the highest element n is inserted into vr', in precisely the position in 
which we inserted it (and we have covered all possible positions in front of elements in the 
active family of tt'), then the result is the primitive permutation tt, as we described it in the 
algorithm. In the first two cases, the active family n enters the active family (and the paired 
value is positive and the Dyck path is rising). In the last two cases, n does not fit into the 
active family, so it removes its friend b respectively r, to start a new family as a bud. There 
are "choice choices of friends, and "choice is recorded as a decoration on the descending Dyck 
path, while in inds, the last entry —k is paired with the "choice choices among the remaining 
unpaired positive entries in inds preceding it. These last observations show that the last entry 
of inds', can be recovered from the pair vr' and tt. As nReds is the number of active ascent 
values, or reds, in vr, this shows that the algorithm is reversible, i.e., the paired values, or 
decorated Dyck path, together with nReds, can be recovered out of vr. Thus the algorithm is 
invertible. 

This ends the proof that primitive permutations are in bijective correspondence with 
paired values, or decorated Dyck paths, together with a specified number of reds and blues in 
the active family. The algorithm also gives an fast and explicit construction for this bijection. 
The smallest cases of data and primitive permutations are the following. 



"Reds 


"Blues 


"Buds 


Paired values 


Dyck paths 


Primitive permutations 


1 


1 





(1,2) 


(+,+) 


(2 1) 


2 


1 





(1,2,3) 


(+,+,+) 


(2 3 1) 


1 


2 





(1,2,3) 


(+,+,+) 


(3 2 1) 


3 


1 





(1,2,3,4) 


(+,+,+,+) 


(2 3 4 1) 


2 


2 





(1,2,3,4) 


(+,+,+,+) 


(3 4 2 1) 


1 


3 





(1,2,3,4) 


(+,+,+,+) 


(4 3 2 1) 


1 


1 


1 


(1,2,3,-1) 
(1,2,3,-2) 
(1,2,3,-3) 


(+,+,+,-(1)) 
(+,+,+,-(2)) 
(+,+,+,-(3)) 


(3 2 (4 1)) 
((4 2) 3 1) 
(2 (4 3) 1) 



Remark that besides the paired values or Dyck paths, one needs to specify also the num- 
ber of reds (ascent values) and blues (descent values) desired, in order to obtain the primitive 
permutation. The structure of the primitive permutation changes considerably with the spec- 
ification of the colors. 

The number of primitive permutations of each type depends only on the sum npam = 
"Reds + "Blues, rather than on nReds and nReds individually, as can be easily checked in the 
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algorithm above. These numbers, in smah cases, are given in the table below. 

nBuds 








1 


2 


3 


4 


5 


2 




3 


21 


207 


2529 


36243 


3 




7 


69 


843 


12081 


197127 


4 




12 


159 


2388 


40221 


751032 


5 




18 


309 


5628 


110781 


2361222 


6 




25 


540 


11760 


268365 


6495345 



npa 



The first line of the table is the sequence oeis A167872: 1, 3, 21, 207, . . ., which is a se- 
quence of moments connected with Feynman numbers. Its double is the sequence A115974: 
2, 6, 42, 414, . . ., which counts the number of Feynman diagrams at perturbative order n which 
are proper, i.e., cannot be obtained from simpler diagrams by concatenation. Note that our al- 
gorithm provides a direct construction. For Feynman diagrams, the known constructions give 
all diagrams, counted by the sequence A000698: 2, 10, 74, 706, . . .. From these, the improper 
diagrams have to be removed case by case. 

The other lines of the diagram, as well as the table itself, do not appear to have been 
studied before. 

11.2 The generating function for primary permutations 



Let 



^{t) 



1 / 1 



2 \t 1 



Er=o(2fc-l)!!t^- 



1 + 3t + 21t^ + 207t^ + 



which is the generating function for the number of primitive permutations with small families 
('^fam = 2) and 0, 1,2, . . . buds, related to the count of fermionic Feynman diagrams. This 
function counts our pairing data for constructing primitive permutations, as {2k — 1)!! counts 
the number of ways to make pairs out of 2k elemnts. 
Let 

ipix) = e^ dy, 

Jo 

which is up to a normalization the imaginary error function. The generating series 

-X + ( e 2^ - 1 (1 - Mt)) + 



a-^r 



+2e 2r 
with the boundary condition 



v- 



l-x 



ij 



1 \\ 1 -t-2tip{t) 



2t 



limx ^^(x,t) = (p(t) 

x—>-0 
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which satisfies the equation 

ox tx X'^ tx'^ 

This function is obtained as a solution to the differential equation which encodes the 
recurrence relation for families which we described before. 

11.3 The decimal code of a permutation 

Let us start with a permutation vr written in line notation. We associate to it a sequence 
resembling the Dewey code for the books in a library. 

The data for a primitive permutation vr consists of a pairing sequence, e.g 

(+1, +2, +3, -2, +4, -1, +5, +6, -4) 

together with colors like (2, 1), i.e. 2 reds and 1 blue, for the unpaired numbers. 

Let us rewrite the data as follows. The paired elements (+4, —4) become 4.1 and 4.2. The 
ending .2 stands for the negative sign. The unpaired elements are followed by their color r or 
b, with all blues before all reds. 

The previous sequence becomes 

(1.1, 2.1, 3.b, 2.2, 4.1, 5.r, 6.r, 4.2). 

We shall call this a colored pairing sequence. The primitive permutation algorithm pro- 
duces a primitive permutation, in which the values, from lowest to highest, correspond exactly 
to our data sequence (1.1, 2.1, 3.b, 2.2, 4.1, 5.r, 6.r, 4.2), in increasing order of the values. 
The colors r and b become the colors indicating ascent values and descent values. 

A bud pair, for instance (+4,-4), now (4.1, 4.2), is in the tree structure of the permu- 
tation, the root, the smallest elements of a new branch. Thus (4.1,4.2) become in the new 
branch the nodes (1, 2), part of a new branch node which is described as before by a primitive 
permutation, e.g 

(1.1, 2.b, 1.2, 3.b, 4.r). 

The numbers 1. • • • 2. • • • which form its root, are (1.1, 2.b), correspond to the bud (4.1, 4.2) 
in the previous colored paired sequence. To unify the notation (4.1,4.2) with (1.1, 2. b), we 
call both of them (4.1.1, 4.2.b). 

We give the prefix 4 to the whole new branch and write it as 

(4.1.1, 4.2. b, 4.1.2, 4.3. b, 4.4.r). 

The previous sequence is now modified to 

(1.1, 2.1, 3.b, 2.2, 4.1.1, 5.r, 6.r, 4.2. b). 

Finally, we prepend to each element the number of its block. Recall that the blocks of the 
permutation are, in the line notation, the segments between successive minima. 
Mark a singlet, which is the block number i, by i.k. 
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Remark Singlets are singlet blocks. At this point, note that each decimal number which 
we have obtained corresponds throughout the process to a value of the permutation. They 
are either first or preceded by the minimum of the previous block, which is smaller. So they 
are always ascent values, to be colored red. By the inverse cycle transform, they correspond 
to on diagonal elements, which are colored black. So depending on the application, singlets, 
or singlet blocks, may be all colored red or black i.k (K is the label of black for computers). 

This way, every value of the permutation gets a unique decimal number separated by 
periods, with all elements numbers except for the last element which is a color. 

In the permutation, a decimal number will correspond to a value of the permutation, 
equal to the position of that decimal number in the list obtained from the permutation. The 
decimal numbers of a primitive permutation will appear in decreasing order. We therefore 
order all decimal numbers in the order of the corresponding values of the permutation. 

The resulting sequence will be called the decimal code of the permutation. 

We shall show that it encodes the permutation we started from, and we shall axiomatize 
it in what follows. Remark at this point that the decimal code has a structure which is very 
far from the one of a permutation. As a permutation, it is closest to vr^^ o vr = id, since 
we used the values of the permutation in ascending order 1,2,3, . . .. Most of the features 
of the decimal code resemble the identity permutation. The subsequence corresponding to a 
primitive node, for instance, has most numbers in natural ascending natural order, except for 
a few pairings to previous numbers. The colors in such a branch are also in ascending order 
from blue to red. The blocks of the permutation which give the first digit of the code are in 
ascending order as well. 

All of the information which allows us to recover the permutation comes from the pairings 
between elements • • • .!.••• .2. • • • , as well as from the shuffling between branches which have 
different heads. Thus, the decimal code is a translation between the order structure of a 
permutation and a pairing and shuffle structure. The latter is typical in many different areas, 
starting from particle physics, with creation paired with annihilation, and time shuffles of 
particles, to biological growth and to the very growth of libraries by addition and loss of 
books, which inspired the Dewey code. 

We shall axiomatize the decimal code structure as follows. 

Definition 

A decimal number i^ is a sequence 11.12- ■ ■ ■ -ik-c separated by dots, in which ii, . . . , ik G 
{1,2,...} and the last entry c is a color, red r or blue b. Any cutoff ii.i2- ■ ■ ■ -ii to length 
/ is called a head of the decimal number i^. A numerical sequence ii.i2---- -ii where 
ii, . . . ,il £ {1, 2, . . . } is called a decimal head. 

A decimal code J is a finite sequence of decimal numbers with the following properties. 

Any decimal head h of an entry of 6 is immediately followed, among the entries of 6, 

(i) either by a color c. In that case, h appears only once as head, h.c is called a leaf, and 
h is called a leaf head. 

(ii) or by a number. In that case h is called a branch head, and is always immediately 
followed only by numbers in the entries of 6. 

Define the node cut of the branch head h as the entries obtained by 

(i) selecting the entries of 5 with head h, in order, 

(ii) cutting off each such entry to length 2 + the length of h, and 

(iii) keeping only the ones ending in 1, 2, or a color, in which each h.n.n' is continued as 
h.n.n' .1. ■ ■ ■ .l.c with a number > of I's after n' . 
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Cut now off the head h and keep the following two digits from each number, in the order 
of the decimal code 5. We ask these numbers to form a colored pairing sequence, i.e., such 
that (i) the numbers of the form n.l, n.b, n.r, should appear in the natural order of n, i.e., 
one of each with n = 1, n = 2, n = 3, . . . in ascending order, 

(ii) each element of the form n.2 should be paired with an element n.l which appeared 
before it, with at most one such pair for each n, 

(iii) the colors c in the elements n.c, which we call leaves, should appear in ascending 
order, first blues then reds, with one of each color, and 

(iv) the sequence of signs obtained by replacing the elements of type n.2 by —1 and n.l, 
n.b, n.r by +1 should form a sequence with the property that the sum of the first k terms is 
> 2 for any A; > 2 (a form of a Dyck sequence). In addition, we require the following property. 

The first digits in each decimal numbers in the decimal code should form a set {1,2,..., m}. 
We call these numbers the block labels. 

We require that for each block label i, there should exist exactly one element a{i), called 
the block anchor, of the form a{i) = i.l. ■ ■ ■ .l.c, with c a color and with a number > of I's 
after i. The color c is black if there are no I's after i, and we call this kind of a{i) a singlet. 
Otherwise, the color c is red. We require the block anchors to appear in ascending order in 
the code 6, so 

<5 = (---,a(l),---,a(2),--- ,a(3),---). 

From the previous conditions a(l) will be on the first place of 6. 
This ends the axiomatization of the code. 

Theorem 11.1. Decimal codes consisting of a sequence of n decimal numbers are in natural 
bijection with permutations of length n. 

The entries of the decimal code are in bijection with the permutation vr correspond to the 
values of IT in ascending order, i.e. to the line form of the inverse permutation vr"^. 

The decimal number corresponding to a permutation value indicates its position in the 
family tree, followed by its color (red or blue for ascent or descent.) The elements with a 
same branch head h correspond, in increasing order, to the values of the elements enclosed 
in a parenthesis in the parenthesized form of vr, which form a branch of the tree form of vr. 
The elements in the node cut of a branch head h form the pairing sequence of the primitive 
permutation corresponding to the root node of the branch of h. The elements with a head i of 
length 1 are the values in the i-th block of tt, the i-th tree in the forest of the tree form of tt. 

Proof. Given a decimal number, which always ends in a color, we call a head of it the decimal 
number obtained by cutting off a number of elements at its end. Cutting off the last entry, the 
color, leaves heads which will be called leaf heads. A leaf head is continued in the decimal 
code by only one color. Cutting off two or more elements at the end leaves nonempty heads 
which will be called branch heads. Branch heads correspond to primitive permutations, 
and to the parentheses in the parenthesized form of the permutation. The algorithm for 
constructing a permutation vr from a decimal code proceeds as follows. 
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For a branch head h, find aU its continuations in the given decimal code which have length 
> \h\ + 2, where \h\ denotes the length of h. These continuations will be of the form h.n.c, 
with c a color, which we keep, and h.n,n', where n,n' are numbers. We keep only those for 
which n' = 1 or 2 and in which the corresponding full decimal number is h.n.n'.l. ■ ■ ■ .l.c, 
where c is a color and the number of I's after h.n.n' is > 0. We now cut off the head h 
and leave the pairs n.c or n.n' with n' = 1 or 2, in the order in which the corresponding full 
decimal numbers appear in the decimal code 6. Our main requirement for a decimal code was 
that such a sequence of pairs is a colored pairing data which encodes a primitive permutation 
o", which we construct with the primitive permutation algorithm. The outer parentheses for 
the permutation a are labeled by the head h. The elements of the active family in a are 
labeled by decimal numbers of the form h.n.c, where c is a color encoding the ascent/descent 
of that member of the active family. The buds of a correspond to paired elements in the data 
sequence and their entries are labeled by pairs h.n.l and h.n.2. Now, assuming that h has 
length > 2, let h = h' .p, where p is a number. The primitive permutation a' corresponding 
to h' will have in it a bud labeled {h' .p.l,h' .p.2) = {h.l,h.2). We now replace this small 
parenthesis in a' with the whole parenthesis containing a, which as we recall was labeled by 
h and contained as a root precisely the pair (h.l, h.2). 

The algorithm proceeds as follows. 

The block anchors are the elements a{i) = i.l. ■ ■ ■ .l.c, where c is a color and the number 
of I's from the second place on is > 0. We cut them off to their first two digits, which will be 
called a{i) and place them in parentheses, to form a list 

((a(l)),...,(a(m))), 

where m is the number of blocks. We call this list the parenthesized decimal code. 

Take now all branch heads of length 1, which label the non singlet blocks. Their two-digit 
extensions, as before, encode primitive permutations. Replace in the parenthesized decimal 
code the singlet (a(i)) with the primitive permutation corresponding to the head i. Continue 
then with branch heads of length 2 and bigger, inductively by head length. 

Construct the primitive permutation corresponding to each branch node head h. Replace 
now in the parenthesized decimal code the bud labeled {h.l, h.2) by the whole parenthesis of 
the primitive permutation built from the head h. We continue to enlarge the parenthesized 
decimal code until all the branch heads have been used and all the corresponding primitive 
permutations have been introduced in it. We now replace each entry in parenthesized decimal 
code, which is a decimal number, by the number indicating the position of that decimal 
number in the list which is the given decimal code 6. We obtain this way the permutation 
vr corresponding to 5, in parenthesized form. Removing the parentheses leaves vr in line 
notation. The parentheses correspond to the nodes which are primitive permutation of the 
tree structure of vr, and are precisely the primitive permutations encoded in the decimal code 
5. The colors at the end of each decimal number give the ascent/descent/singlet statute of 
the corresponding permutation value. It is not hard to see that the two directions of the 
algorithms matching a permutation with a decimal code are inverse to each other. The main 
difficulties are concentrated in the primitive permutation algorithm, which we described in 
detail. D 
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An example The code for the permutation 

(3, 16, 9, 19 ,1, 20, 4, 18, 10, 13, 7, 11, 14, 8, 12, 6, 2, 15, 17, 5) 
in which the colors indicate ascent /descent values, is 



1 


l.l.b 


2 


2.1.b 


3 


1.2.r 


4 


2.2.b 


5 


3.1.b 


6 


2.3.1. Lb 


7 


2.3.2.1.b 


8 


2.3.3.b 


9 


1.3.1.b 


10 


2.3.2.2.1.b 


11 


2.3.4.r 


12 


2.3.1.2.r 


13 


2. 3. 2.3. r 


14 


2.3.5.r 


15 


3.2.r 


16 


1.3.2.r 


17 


3.3.r 


18 


2.3.2.2.2.r 


19 


1.4.r 


20 


2.4.r 



The block anchors are a(l) = l.l.b, a(2) = 2.1.b, a(3) = 3.1.b. The parenthesized 
decimal code starts with the two digit heads of these as 

((1.1), (2.1), (3.1)). 

The branch heads, which correspond to primitive permutations are obtain by cutting away 
> 2 digits at the end of the decimal numbers. They are the following. 

1 
2 
3 

1.3 

2.3 

2.3.1 

2.3.2 

2.3.2.2 

The decimal numbers on the branch with head h = 2.3, selected and truncated to the two 
digits following h, are in order 

1.1 2.1 3.b 2.2 4.r 1.2 5.r. 
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In pairing notation, they correspond to (1,2,3,-2,4,-1,5). The corresponding primitive 
permutation is ((4, 2), 5, 7, 3, (6, 1)), with active family (5, 7, 3). The buds are (4, 2) and (6, 1). 
Recalhng that the decimal head was h = 2.3, the first four digits of the decimal code of the 
numbers in each parenthesis are 

((2.3.2.2, 2.3.2.1), 2.3.4.r, 2.3.5.r, 2.3.3.b, (2.3.1.2, 2.3.1.1)). 

Note that the head h = 2.3 encodes this whole primitive permutation, and in particular 
labels its outer parentheses. The inner buds, and their parentheses are encoded by 2.3.2 and 
2.3.1, respectively. The endings (.2, .1) in each indicate the members of the bud, which will 
become the root elements, cofounder and founder, of the branches in which these buds would 
develop. These emerging branches have heads 2.3.2 and respectively 2.3.1. The branches will 
be inserted in the parenthesized decimal code, in the place of these buds. 

The last two digits of these, which encode the values of the primitive permutation 
((4, 2), 5, 7, 3, (6,1)) are 

((2.2,2.1),4.r,5.r,3.b,(1.2,l.l)). 

In the end, we obtain a parenthesized decimal code which contains all the entries in the 
given decimal code and starts with 

((1.2.r, 1.3.2.r, 1.3.1.b, 1.4.r, •••))• 

We replace these numbers by their position in the decimal code 5, to obtain the parenthe- 
sized 

((3, (16, 9), 19 ,1), (20, 4, (((18, 10), 13, 7), 11, 14, 8, (12, 6)), 2), (15, 17, 5)). 

This is precisely our initial permutation. The branch with head h = 2.3 which we analyzed 
before, called there ((4, 2), 5, 7, 3, (6, 1)), have now become 

(((18,10), 13, 7), 11, 14, 8,(12, 6)), 

in which the bud (4, 2) has expanded to the branch ((18, 10), 13, 7). We keep out of this 
branch its smallest elements (10, 7), which form its root as cofounder and founder. That 
branch started in the permutation when 10 pulled out 7 from its family to start a new branch. 
Replacing that branch by its root, we obtain the primitive node 

((10, 7), 11, 14, 8, (12, 6)). 

The pattern of this node, obtained by reducing the numbers, in order, to the smallest 
possible values, is 

((4, 2), 5, 7, 3, (6, 1)), 

which is the primitive permutation of the node. This primitive permutation is inserted into 
the final parenthesized permutation by the presence of the head h = 2.3 which labeled it in 
the decimal code. 

The decimal code is shown this way to encode the structure of the permutation vr by means 
of the values corresponding to the pairings in each primitive component, together with the way 
in which the values corresponding to different branches are shufSed with each other. Other 
than this, the components of the decimal code are all in natural ascending order. The order 
change described by a permutation has been replaced with pairing and shuffling information, 
of an entirely different nature. 
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12 Family statistics: modeling cell biology 

Just like the combinatorics of binomial coefficients, the 2" subsets of n into subsets with the 
same cardinality, led to the normal distribution and the start of modern probability theory, we 
expect the map from of each of the n! permutations of n into trees of primitive permutations 
and into registries of families to open new directions in probability theory. The fact that 
families only grow when they are very small, but overwhelmingly bud into new branches when 
they are large, and the fact that each branch grows independently, with branches "competing" 
for new numbers, could be ideal for modeling biological processes. 

The n elements of a line permutation are partitioned into families playing the role of cells 
with k blue and / red elements each, counting the ascents and descents. As new numbers, 
modeling food, are introduced into the permutation, cells receive it proportional to their size. 

Out of the k + I positions of a {k,l) cell, l/{k + I) (or 2/{k + I) in case k is equal to 
1), produced cell growth. The other positions diminish the cell by one and produce a (1, 1) 
offspring. Thus in early stages cells mostly grow and when they become larger they mostly 
reproduce. 

We have computed the statistical distribution of cells by (k, I) type and the result was 
remarkably simple, a multinomial with multivariate Gaussian distribution. 

The number of permutations of n with ni singlets and ?T-(i j) families with i ascents and j 
descents is 

n\ (E"-(i,j))! 

The configurations are all the sequences of numbers (ni , f^(i,i) , ^{2,1) ; '^{1,2) ; '^(3,1) > "-(2,2) > ^(i,3) > • • 
with 

ni + 2n(^i) + 3(n(2,i) + n(i_2)) + 4(n(3^i) + re(2,2) + ^(1,3)) + ■■■ =n 

Summed over these the multinomials above add up to n\. The reason behind this very simple 
form of the sum at the end of a very complex process is that our structure theorem gives a 
bijection between permutations of n and unrestricted compositions of n into families. 

Thus a complex cell growth model has its outcome exactly described, in a simple way, 
by our new structure theory. It opens a gate towards the study of the deeper structure of 
the model, the possible tree structures arising from it, and corresponding new structures in 
probability and statistics. 

13 Conclusions 

In this paper, we have developed aspects of the internal structure of permutations, as well as 
the algebraic structure of multinomials, which count permutations by means of bicolored set 
partitions. 

In shortly forthcoming papers, we shall describe applications of these results in simplicial 
geometry and in theoretical physics. 
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